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ABSTRACT 

The aim of this paper is to set constraints of the epochs of early-type galaxy formation through the 
'archaeology' of the stellar populations in local galaxies. Using our models of absorption line indices 
that account for variable abundance ratios, we derive ages, total metallicities, and element ratios of 124 
early-type galaxies in high and low density environments. The data are analyzed by comparison with 
mock galaxy samples created through Monte Carlo simulations taking the typical average observational 
errors into account, in order to eliminate artifacts caused by correlated errors. We find that all three 
parameters age, metallicity, and a/Fe ratio are correlated with velocity dispersion. We show that 
these results are robust against recent revisions of the local abundance pattern at high metallicities. 
To recover the observed scatter we need to assume an intrinsic scatter of about 20 per cent in age, 
0.08 dex in [Z/H] and 0.05 dex in [a/Fe]. AU low-mass objects with < 10^° Mq (ct < 130 km/s) 
show evidence for the presence of intermediate-age stellar populations with low a/Fe ratios. About 
20 per cent of the intermediate-mass objects with 10^° < M^/Mq < lO^i (110 < cr/km/s < 230, both 
ellipticals and lenticulars) must have either a young subpopulation or a blue horizontal branch. Based 
on the above relationships valid for the bulk of the sample, we show that the Mg-cr relation is mainly 
driven by metallicity, with similar shares from the a/Fe ratio (23 per cent) and age (17 per cent). 
We further find evidence for an influence of the environment on the stellar population properties. 
Massive early-type galaxies in low-density environments appear on average ^ 2 Gyrs younger and 
slightly (^ 0.05 — 0.1 dex) more metal-rich than their counterparts in high density environments. No 
ofl'sets in the a/Fe ratios, instead, are detected. With the aid of a simple chemical evolution model, 
we translate the derived ages and a/Fe ratios into star formation histories. We show that most star 
formation activity in early-type galaxies is expected to have happened between redshifts ~ 3 and 5 in 
high density and between redshifts 1 and 2 in low density environments. We conclude that at least 50 
per cent of the total stellar mass density must have already formed at z ~ 1, in good agreement with 
observational estimates of the total stellar mass density as a function of redshift. Our results suggest 
that significant mass growth in the early-type galaxy population below z ~ 1 must be restricted to 
less massive objects, and a significant increase of the stellar mass density between redshifts 1 and 2 
should be present caused mainly by the field galaxy population. The results of this paper further 
imply vigorous star formation episodes in massive objects at z ~ 2-5 and the presence of evolved 
ellipticals around z 1, both observationally identified as SCUBA galaxies and EROs, respectively. 
Subject headings: galaxies: abundances - galaxies: elliptical and lenticular, cD - galaxies: stellar 
content ~ galaxies: formation - galaxies: evolution 



1. INTRODUCTION 



The most direct way to constrain the formation and 
evolution of galaxies certainly i s to trace back their evo- 
lution with redshift (e.g. lAragon-Salamanca et al.l 
19931 iB ender et al.l 119961; iKelson et al.' I200C " 
K odama e t al. 19 99; J oggianti et al. 1999; Sagl ia et al.l 
2p00; Stanfor d efairi l998: van Dokkum et ain2f 
Ziegler fc Bended ll997t .Ziegler et al . 1999). The price 
to be paid, however, is that high-redshift data naturally 
have lower quality and are therefore more difficult 
to interpret. A clear complication is the so-called 
progenitor bias, which implies that galaxies observed 
at low and high reds hift are not necessarily drawn 
from the same sample ijvan Dokkum et ahl 120001) . The 
alternative approach is the detailed investigation of 



the stellar populations in local galaxies, which has 
been pioneered by analyzing slopes and scatter of 
color-magnitud e and scaling rela t ions o f early-type 
galax i es (e.g. 'Bender ct al. 1992", '199?; 'Bower et aU 
19921 iDiorgovski fc Davis, .1987; .Dressier et al.i 119811 
Renzini fc Ciottil Il993() . followed by a number of 
detailed studies of absorption line indices (e.g 



Pglgtisil 119891 iGorgas. E fstathiou. fc Aragon Salaraaj 



m 



iWorthev. Faber. fc Gonzaled 



ijavies. Sadler. fc Peletieil ll993t iCarollo fc Danzigei 
19941 iRose et all 119941 IBender fc Paauet 



iFisher. Franx. fc lUin gworthf 

Jorge ns en. Franx . fc Kia^rgaard 199 5t [ Greggiol 

yazdekis et al.l 119 97: Tantalo. Chiosi. fc Br^essar 
iKuntschner fc Davies, .1998; .Mehlert et al, 
19981 iWorthevI 119981 iKobavashi fc Arimotc 
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2003al and others). We call this the 'archaeology 
approach'. The confrontation with predictions from 
models of galaxy formation is certainly most mean- 
ingful, when the two approaches, the mining of the 
high-redshift universe and the archaeology of local 
galaxies set consistent constraints. In this paper we 
follow the latter approach. 

The main challenge in the archaeology of stellar pop- 
ulation is the disentanglement of age and metallic- 
ity effects. The use of ab sorption line indices to lift 
this degeneracy is pow erful (|Faber et alJll985t iGonzale^ 
119931 lWorthevlll994() . but has been up to now ham- 
pered by the fact that different metallic line indices 
yield different metallicities and therefore different ages 
(see references above). We have solved this prob- 
lem by developing stellar population models that in- 
clude element abundance ratio effects and now allow 
for an un-ambiguous derivation of age, total metallicity, 
and element ratios from (Lick) absor ption line indices 
ijThomas. Maraston. fc Bendeill2003bl hereafter TMB). 
A complication that remains, however, is the degeneracy 
between age and horizontal branch morphology, which 
stems from the fact that the presence of warm horizontal 
branch stars (which cannot be excluded) strengthens the 
Balmer ab sorption and can m imic a younger st ellar popu- 

fc Barbuv 199,4 lLee et all 
We will additionally 
discuss this problem in the present paper. 

The principal aim of this study is to constrain the for- 
mation epochs of the stellar populations in early-type 
galaxies as a function of their type (elliptical and lentic- 
ular), mass and environmental density. For this purpose 
we analyze a homogeneous, high-quality data sample of 
124 early- type galaxies in various environmental densi- 
ties. With our new stellar population models (TMB) we 
determine ages, total metallicities, and a/Fe ratios from 
the absorption line indices H/3, Mg6, and (Fe), and seek 
for possible correlations of these parameters with veloc- 
ity dispersio n. First tentative results have alread y been 
discussed in iThomas. Maraston. fc Bendeil l)2002i) . Here 
we present the final and more comprehensive analysis. In 
particular, we compare the observed data set with mock 
samples produced through Monte Carlo simulations tak- 
ing the typical average observational errors into account, 
in ord er to eliminate the confusion caused by correlated 
errors llKuntschner et 111120011: ITerlevich fc Forbe^l2002l: 
iTrager et al.ll2000aD . 

The a/Fe element ratio plays a key role for the ac- 
complishment of our main goal, namely the derivation of 
formation epochs. While the so-called a-elements O, Ne, 
Mg, Si, S, Ar, Ca, Ti (particles that are build up with 
a-particle nuclei) plus the elements N and Na are deliv- 
ered mainly by Type II supernova explosions of massive 
progenitor stars, a substantial fraction of the Fe-peak 
elements Fe and Cr co mes from the delayed exploding 
Type la supernovae fe.g. lNomoto. Thielemann, fc Yokoil 



19841: iThielemann. Nomoto. fc Hashimotol 119961: 
Wooslev fc Weaved Il995j) . Hence, the a/Fe ra- 
tio quantifies the relative importance of Ty pe II 
and Type la supernovae iGreggio fc Renzinil 
Matteucci & Grcggio 1986'; Pag el fc Taiitvaisienel 



1983; 



Thomas , _Grej gio. fc Bender 199^ . and therefore carries 



information about the timescale over which star forma- 
tion occurs. Thus, the a/Fe ratio can be considered as 
an additional measure of late star formation, and we 
will use it both to constrain formation timescales and to 
lift the degeneracy between age and horizontal branch 
morphology. 

The paper is organized as follows. In Sections El and |31 
we present the methodology of the paper, i.e. we intro- 
duce the data sample, we briefly summarize the stellar 
population model of TMB, and we explain the construc- 
tion of the mock galaxy sample through Monte Carlo 
simulations. The results of this analysis are presented in 
Section^ We use these results to explore the origin of 
the Mg-(T relation (Sectional and to derive star forma- 
tion histories (Section with the aim to constrain the 
epoch of early-type galaxy formation. The main results 
of the paper are discussed and summarized in Sections [3 
and El 

2. THE DATA SAMPLE 

2.1. Sample selection 

We analyze a sample of 124 early-type galaxies, 70 of 
which reside in low-density and 54 in high-density en- 
vironments, containing roughly equal fractions of ellip- 
tical and lenticular (SO) galaxies. To maintain statisti- 
cal significance, in the following we make the rough bi- 
modal distinction between low-density and high-density 
environment. It should be emphasized that the clas- 
sification 'low density' comprises galaxies in environ- 
ments with low density, not necessarily isolated objects. 
The sample is constructed from the following sources: 
41 galaxies i n high (Virgo c luster) and low density en- 
vironments fGonzalez"1993), 32 Coma cluster galaxies 
(Mchlcrt ct al. 2000, 2003), a nd 51 galax i es ( mostly low 
density environment) from Bcuing et al. (2002) selected 
from the ESO-LV catalog (^Lauberts fc Valentiin .1989.) 
(highest-quality objects). In this latter sample, objects 
with a local galaxy surface density NGt > 9 are as- 
signed to the high density en vironment. NGt is given in 
iLauberts fc ValentiinI (|1989) and is the number of galax- 
ies per square degree insid e a radius of one deg ree around 
the considered galaxy fsee lBeuing et al.ll2002fl . The sam- 
ple spans a large range in central velocity dispersion 
60 < aa /kms~^ < 340 with similar distributions in both 
environments (Fig.^. 

In the data samples quoted above, absorption line 
strengths are measured as functions of galaxy radius. We 
adopt the central indices determined within 1/10 of the 
effective radius, so that the analysis presented here does 
not suffer from aperture effects. We use the Lick indices 
H/3, Mgb, and (Fe) = 0.5(Fe5270 Fe5335) to derive 
the stellar population parameters age t, total metallicity 
[Z/H], and the element abundance ratio [a/Fe]. The me- 
dians of the l-cr errors in H/3, Mg6, Fe5270, and Fe5335 
are 0.06, 0.06, 0.07, 0.08 A, respectively. 

We have re-observed 19 objects of the iBeuing et alJ 
1)2002(1 sample with the same telescope (the 1.5m tele- 
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Fig. 1. — Distribution of central velocity dispersion in high den- 
sity (grey shaded histogram) and low density (hatched histogram) 
environments. 



scope on La Silla) with the principal aim of extending 
the wavelength range to the blue (Mendes de Oliveira 
et al., in preparation). We find a good agreement for 
the optical absorption line indices H/3, Mg6, and (Fe) 
between the two samples, with a small systematic offset 
in Mg6 of 0.2 A. As the latter observations have higher 
quality due to the installation of a new CCD chip and 
holographic grating, we c orrected for the offse t in Mg5 
and replaced the values of iBeuing et alJ f2002) with the 
new measurements for the objects in common. We em- 
phasize here, however, that the results of this paper are 
not affected by this procedure. 

2.2. Absorption line indices 

In Fig. |21we present the absorption line index measure- 
ments plotting [MgFe]' vs. H/3 (top panels) and Mg6 vs. 
(Fe) (bottom panels). The index 

[MgFe]' = v'Mg b ■ (0.72 • Fe5270 -I- 0.28 • Fe5335) , 

a slig ht modification of [MgFe] defined by IGonzaled 
()1993j) . is almost completely independent of a/Fe ratio 
variations (TMB). Left and right panels show data in 
high density and low density environments, respectively. 
Red symbols are elliptical, blue symbols are lenticular, 
and green symbols are cD galaxies. Overplotted are our 
stellar population models (TMB) for various ages, metal- 
hcities, and a/Fe ratios as indicated by the labels in the 
diagrams. 

2.2.1. TMB models 

The TMB models take the effects from element abun- 
dance ratio changes on all Lick indices into account, 
hence give Lick indices of simple stellar populations 
(SSPs) not only as functions of age and metallicity 
but also as a function of the a/Fe ratio. They are 
based on t he ev o lution ary population synthesis co de of 
iMarastod ifTflM l200l . In iThomas eTaP l|2003bD . the 
impact from ele ment r atio changes is computed with 
the help of the iTripTc co fc Bell (1995) response func- 
tions, using an ext ension of the method introduced by 
iTrager et all l|2000b.) . The updated version used in this 
paper adopts the new metallicity-dependent response 



functions from iKorn. Maraston. fc Thomad l)2005j) . As 
discussed there, the impact of metallicity on the response 
func tions for the indices c onsidered here is negligible (see 
alsolril omas. Maraston. fc Korn 2004), hence does not 
affect the results of this paper. 

Particular care was taken to calibrate the SSP mod- 
els with globular cluster data ijMaraston et alJ [2003) , 
which — most importantly — include objects with rela- 
tively high, namely solar, metalhcities (jPuzia et alJ 
I2002|) . We match very well their H/3, Mg6, and (Fe) 
indices with models in which the a/Fe ratio is enhanced 
by a factor two relative to the solar value, in agreement 
with results from high-r esolution spectroscopy of individ- 
ual stars (see TMB and lMaraston et alJl200ti and refer- 
ences therein). Because of the inclusion of element ratio 
effects, the models allow for the clear distinction between 
total metallicity [Z/H] and the a to iron-peak elements 
ratio [a/Fe]. 

We use these models to derive the three stellar popu- 
lation parameters age, metallicity, and a/Fe ratio from 
the three line indices H/3, Mg&, and (Fe) in a twofold 
iterative procedure explained below. As mentioned in 
the Introduction, there is a degeneracy between age and 
horizontal branch morphology. A way to lift this degen- 
eracy is the consideration of the a/Fe ratio. We expect 
young stellar populations generally to be characterized 
by lower a/Fe ratios owing to the late enrichment of Fe 
from Type la supernovae. 

In this paper we present a detailed investigation of the 
data set based on Monte Carlo simulations as shall be 
discussed in the following sections. However, first tenta- 
tive conclusions can already be drawn by the direct com- 
parison between the observational data and SSP model 
grids, as discussed in the following. 

2.2.2. Ages and metallicities 

The index [MgFe]' is a good tracer of total metallic- 
ity (still degenerate with age, of course), as it is almost 
completely independent of a/Fe ratio variations (TMB). 
As also H/3 is only very little sensitive to a/Fe (TMB), 
the top panels of Fig. |21 are well suited to read off ages 
and total metallicities independent of abundance ratio 
effects. 

The sample splits in two subclasses divided by the 
Balmer absorption index at H/? « 2 A. The main part 
of the data is centered around H/3 — 1.6 A, while 
about one quarter of the objects have relatively high 
H/? absorption implying the presence of either young 
stell ar populatio ns or blue horizontal branch morpholo- 
gies (jMaraston fc Thomas 2000). This su rprisingly well- 
define d separation is also discussed in IMehle rt et alJ 
ll2003ll for the Coma cluster objects alone. Interestingly, 
the pattern is still present for the larger galaxy sample 
considered here, and in particular the fraction of high H/? 
objects seems to be independent of the environmental 
density (Fig. but we do not find a significant differ- 
ence with respect to the galaxy type. Fig. |21 further hints 
toward a difference between objects as a function of envi- 
ronment within the old (low H/3) class, in the sense that 
galaxies in low density environments seem to have on av- 
erage slightly younger ages and higher metallicities. This 
difference is more clearly visible in the direct comparison 
of the derived ages discussed below. 
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Fig. 2. — Lick indices [MgFe]' versus H/3 ( top panels) and Mg b vers us (Fe) (bottom panels). Symbols are the early- type galaxies analyzed 
in this paper, taken from Gonz aIe3 119931) . iBeuing et al.l 1120021 ). and IMehlert et ahl 120031) . Index values are measured within ^ 1/10 Tg. 
Left and right panels show objects in high and low density environments, respectively. Red symbols arc elliptical, blue symbols lenticular, 
and green symbols are cD galaxies. Median 1— ct error bars are shown. SSP models (TMB) with the metallicities [Z/H] = 0.0, 0.35, 0.67 are 
plotted for different ages t = 2, 3, 5, 10, 15 Gyr at fixed a/Fe ratio ([a/Fe] = 0) (top panels) and different a/Fe ratios [a/Fe] = 0.0, 0.3, 0.5 
at fixed age {t = 12 Gyr) (bottom panels) as indicated by the labels. 



2.2.3. Element abundance ratios 

The bottom panel of Fig. [21 can be used to get rough 
estimates for the a/Fe ratios of the sample galaxies. It 
should be emphasized, however, that the Mg6-(Fe) plane 
is not completely degenerate in age, and the ages derived 
from the top panels have to be taken into account for the 
precise derivation of a/Fe ratios. Still we can infer from 
the bottom panels of Fig. El that the a/Fe ratios range 
roughly between solar and two times solar, and that ob- 
jects with strong Mg b indices, hen ce velocity dispersio n 
because of the Mg-cr relation (e.g. IBender et aITll993j) . 
tend to have higher a/Fe ratios. At least no obvious im- 
pact from environmental density on the a/Fe ratios is 
visible. 

2.2.4. Derivation of stellar parameters 

The stellar population parameters are derived in a two- 
fold iterative procedure. First, we arbitrarily fix the a/Fe 
ratio, and determine ages and metallicities for the index 
pairs (H/3,Mg6) and (H/3,(Fe)), by starting with arbi- 
trary age-metallicity pairs, which we modify iteratively 
until both index pairs are reproduced. The two metallic- 
ities obtained from Mg5 and (Fe), respectively, are used 
to adjust the a/Fe ratio, and to start a new iteration. 
These steps are repeated until the age-metallicity pairs 



derived from (H/3,Mg6) and (H/3,(Fe)) at a given a/Fe 
ratio are consistent within 1 per cent accuracy. For ages 
and metallicities between the grid points quoted above, 
we interpolate linearly. 

The resulting ages, metallicities, and a/ Fe fo r the sam- 
ple under investigation are listed in Table IXSl 

3. THE MONTE CARLO APPROACH 

An uncomfortable effect of the degeneracy between 
age and metallicity is that the errors of these stel- 
lar population parameters are not independent. It 
is still a matter of debate, for instance, whether 
the age-metallicity anti-correlation found for early-type 
galaxies in a numbe r of studies is an artifact caused 
by correlated errors iJorgensenlllOO^ iKuntschner et alJ 
20011 i Poggianti et al.' '2001 bt iProctor fc SansomI 120021: 
Terlevi ch & Forbes 2002; Tr ager et al.ll2000a|) . A pos- 
sibility to get a handle on these effects is the perfor- 
mance of Monte Carlo simulations taking into account 
the observational uncertainties of the line indices from 
which the stellar population parameters are derived. 
Kuntschne r et all |2001) follow this method in order to 
assess the relationship between ages and metallicities of 
72 early-type galaxies by comparison with a mock sample 
produced through Monte Carlo simulations. 
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Fig. 3. — Stellar population parameters of an object with 
H/3 = 1.59 A, Mgb = 4.73 A, and (Fe) = 2.84 Acorresponding to 
an age of 10.7 Gyr, a metallicity [Z/R] = 0.26, and [a/Fc] = 0.25 
as indicated by the dotted hnes. The circles are Monte Carlo re- 
alizations perturbing the original index values with the l-cr er- 
rors dH/3 = 0.06 A, dMgb = 0.06 A, dFe5270 = 0.07 A, and 
dFe5335 = 0.08 A). The resulting 1-a errors of the stellar popula- 
tion parameters are 1.48 Gyr in age, 0.04 dex in total metallicity, 
and 0.02 dex in a/Fe ratio. 



In this pape r 
iKuntschner et all ll200 



we 



extend the approach of 
to the full parameter space 
including also a/Fe ratios besides age and total metal- 
licity. We produce a sample of artificial galaxies with 
given ages, metallicities, and a/Fe ratios, from which we 
compute with the TMB SSP models the line indices H/3, 
Mg6, and (Fe). By means of Monte Carlo simulations 
we perturb these 'exact' values with the 1-a errors 
quoted above, assuming a Gaussian error distribution. 
An illustrative example for the resulting correlations 
between the errors of the three population parameters is 
given in Fig.|31for an artificial object with H/3 = 1.59 A, 
Mg6 = 4.73 A, and (Fe) = 2.84 A corresponding to 
an age of 10.7 Gyr, a metallicity [Z/Fl] = 0.26, and 
[a/Fe] = 0.25 indicated by the dotted lines. The circles 
are 100 Monte Carlo realizations taking the 1-cr errors 
in the line indices (dH/3 = 0.06 A, dMgb = 0.06 A, 
dFe5270 = 0.07 A, dFe5335 = 0.08 A) into account. 

The distributions of the resulting stellar population pa- 
rameters can be well described by Gaussians. The cor- 
responding 1-a errors are 1.48 Gyr in age, 0.04 dex in 
total metallicity, and 0.02 dex in a/Fe ratio. We confirm 
that errors produce a very tight linear anti-correlation 
between log age and metallicity (top left panel) as al- 
ready mentioned above. Interestingly, we find that also 
the errors of age and a/Fe ratio are anti-correlated (bot- 
tom left panel). This relationship is the result of the 
slightly stronger age dependence of the Mg b index with 
respect to (Fe), caused by the large r contribution from 
turnoff stars to Mg& and Mg2 (^Maraston et al.ll2003D . 
An underestimation of the age leads to the prediction 



of lower Mg indices relative to iron, which results in the 
derivation of a larger a/Fe ratio, and vice versa. 

In the following Monte Carlo simulations will be used 
to reproduce and analyze the distribution of the data 
in the diagrams shown in Fig. |21 The simulations are 
carried out in two steps, low density and high density 
environments are considered separately throughout the 
analysis. 

1) First, we seek correlations between galaxy veloc- 
ity dispersion and the stellar population properties age, 
metallicity, and a/Fe ratio t hat reproduc e the majority 
of the data sample. Indeed, iTraeer e~l . (2000a) have 
already shown that velocity dispersion is the only struc- 
tural parameter that is found to modulate the stellar 
populations of galaxies. As we aim at deriving average 
properties, simple stellar populations of single age and 
single chemical composition are used. We start with the 
velocity dispersion (a) distribution of the data sample 
shown in Fig. ^ Based on arbitrarily chosen relation- 
ships, we assign to every a an age, metallicity, and a/Fe 
ratio, which we transform into H/3, Mg b, and (Fe). These 
indices are then perturbed with the observational error 
(see above) adopting a Gaussian probability distribution. 
In this way we obtain a mock galaxy sample that is clos- 
est to the observational situation in terms of size, velocity 
dispersion distribution, and error perturbation. If neces- 
sary, we add an intrinsic scatter in the stellar population 
parameters age, metallicity, and a/Fe. This step is re- 
peated, until the majority of the observed data points in 
the index-index diagrams [MgFe]' vs. H/3 and Mg& vs. 
(Fe) is reproduced (comparison 'by eye'). The quality of 
the fit is further checked with respect to the zero points, 
slopes, and scatter of the relationships between the stel- 
lar parameters (age, metallicity, a/Fe ratio) and velocity 
dispersion. By means of the index-index diagrams we 
identify those objects that are not reproduced by the 
mock data sample based on the simple relationships in- 
troduced in this step. 

2) In a second step, we therefore perturb the general 
simple scaling laws found in the first step in 2-component 
models in order to match those remaining objects. Like 
before, we adopt the exact a distribution of the outliers 
and assign age, metallicity and a/Fe to every a based on 
the above final relationships to define the base popula- 
tion. Then we add a second component with arbitrary 
stellar population parameters and determine the index 
values of this composite population. These are then per- 
turbed like in the first step and compared to the obser- 
vational data. We note that the contribution in mass of 
the second component is treated as the forth free param- 
eter, hence the solutions will not be unique, as we have 
only three observables as constraints (H/3, Mg6, (Fe)). 
As discussed in detail below, we will use plausibility ar- 
guments trying to maintain astrophysical feasibility of 
the Monte Carlo realizations, in order to constrain the 
relative weights of the subcomponents. 

4. RESULTS 

In Fig. ^ the comparison of our Monte Carlo realiza- 
tions (black symbols) with the observational data (col- 
ored symbols) is shown in the index-index planes intro- 
duced with Fig. 121 By means of the 'age-metallicity' di- 
agram (top panels of Fig.^J we distinguish between four 
different categories, independent of the environment, as 
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and 1.4 Gyr old subcomponents (10 and 20 per cent in mass, respectively) on top of the above relationships, respectively. Diamonds are 
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indicated by the symbol types: 1) V^-light averaged ages 
older than ^ 5 Gyr (circles), 2) ages between 2 and 3 Gyr 
and metallicities [Z/H] < 0.35 (triangles), 3) ages be- 
tween 2 and 3 Gyr and metallicities 0.35 < [Z/R] < 0.67 
(squares), 4) ages between 2 and 5 Gyr and metallicities 
[Z/R] > 0.67 (diamonds). 

Here we provide a short summary how these categories 
are reproduced, detailed descriptions are discussed below 
in Sections 4.1 to 4.4. Category 1 contains intermediate- 
mass and massive galaxies (120 < cr/(km/s) < 320) and 
can be reproduced by linear relationships between the 
stellar population parameters and log cr. For simplic- 
ity we shall call this category old population, as Cate- 
gories 2-4, unlike Category 1, require 2-component mod- 
els with young sub-populations. Category 2 consists of 
low-mass galaxies (50 < (T/(km/s) < 130) and are repro- 
duced assuming a minor (10 per cent in mass) young 
(~ 1.2 Gyr) sub-component. Category 3 are young 
intermediate-mass galaxies (110 < fT/(km/s) < 230) 
with a minor (20 per cent in mass) young (~ 1.4 Gyr) 
component, while Category 4 contains young massive 
galaxies (190 < (T/(km/s) < 360) with a major (70 per 



Table 1. Model parameters 





circles 


triangles 


squares 


diamonds 


(T (km/s) 


120-320 


50-130 


110-230 


190-360 


log A/*/M0 


10-12 


< 10 


10-11 


> 11 


1st comp. 


Eqn.0 


Eqn.0 


Eqn.0 


Eqn.[I] 


2nd comp. 


no 


yes 


yes 


yes 


Mass 




10% 


20% 


70% 


Age (Gyr) 




1.2 


1.4 


2.5 


[Z/H] 




~ 0.2 


~ 0.8 


~ 0.7 


[a/Fe] 




0.06 


~ 0.20 


~ 0.34 



Note. — Stellar masses (row 2) are from Eqn. |^ The last five 
rows refer to the second component. 



cent in mass) young (^ 2.5 Gyr) component. The data 
sample is well characterized by these four object classes 
and well reproduced by the JVIonte Carlo simulations. In 
the following section we will discuss in detail the ir ch ar- 
acteristics. They are summarized in Table ITI and IA2I 

4.1. Old population - circles 
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Fig. 5. — Absorption line indices as a functions of velocity disper- 
sion (measured within 1/10 r^) for the 'old' subpopulation shown 
as circles in Fig. [H (re d: elliptical, blu e: SO, green: cD, tak e n from 
IGonzaled 11293), Bcu ing et al.l <200l) . and IMehlert et al.l ll2003i . 
measured within 1/10 re). Median I-ct error bars are shown. The 
Monte Carlo realizations are shown by the black symbols. The 
latter take the observational errors (see Fig. |3J into account plus 
an intrinsic scatter A[o/Fe] = 0.05 dex, A[Z/H] = 0.08 dex, and 
At = 0.2 ■ t Gyr in high density environments (left panel) and 
A[a/Fe] = 0.07 dex, A[Z/H] = 0.10 dex, and At = 0.25 ■ t Gyr in 
low-density environments (right panel). 



This subpopulation represents the major fraction of 
the present sample, namely 60 and 70 per cent in 
high density and low density environments, respectively 
containing both elliptical and lenticular galaxies. The 
objects display relatively well-defined relationships be- 
tween absorption line indices and central velocity dis- 
persion, as shown in Fig. by the colored symbols. 
Both metallic indices Mgb and (Fe) correlate with ve- 
locity disp ersion, which is well-known for Mgb (e.g., 
iBender et al. 1993 !), and has also been recently es- 
tabhshcd for (Fc) (Bernardi et all 120031 iCaldwell et"an 
12003: Kun tschncr et al. 2001) . The Balmer index 
H/3, instead , clearly tends to de c rease with increasing 
a (see also iBernardi et alJ l2003t ICaldwell et all l200a 
IKuntschner et al.ll200H) . Fig. shows that these rela- 
tionships are essentially independent of the environment, 
except that the (Fe)-f7 relation seems somewhat flatter 
in low-density environments. 



4.1.1. Scaling relations 

The black symbols in Fig. [3 are the Monte Carlo re- 
alizations corresponding to the dominating 'old' popu- 
lation (circles in Fig. 0J. In these simulations we have 
assumed linear correlations of the parameters age, metal- 
licity and a/Fe ratio with log cr. The dependencies that 
fit best all three indices as functions of velocity dispersion 
for objects in high density environments are (quantities 
in brackets are the values derived for the low-density en- 
vironment): 

[a/Fe] = -0.42 (-0.42) + 0.28 (0.28) log a (1) 
[Z/H]==-1.06 (-1.03) -h 0.55 (0.57) loga 
log VGyr = 0.46 (0.17) + 0.238 (0.32) loga 

The simulations take the observational errors (see Fig.|3Jl 
into account plus an intrinsic scatter A [a/Fe] = 0.05 dex, 
A[Z/H] = 0.08 dex, and At = 0.2 ■ t Gyr in high den- 
sity environments and A [a/Fe] — 0.07 dex, A[Z/H] = 
0.10 dex, and At = 0.25 • t Gyr in low-density environ- 
ments. The relationships with velocity dispersion includ- 
ing their scatter of the line indices are well reproduced. 
The relative contributions of the observational errors and 
the intrinsic scatter of the stellar population parameters 
to the overall spread in the data and simulations will be 
discussed in more detail below. 

Before doing so, we provide estimates of the stellar 
masses as a function of velocity dispersion. For this pur- 
pose we transform the ages and metallicities of Eqn. ^ 
into stellar mass-to-light ratios with the help of stel- 
lar population models a nd use these to derive (via the 
Faber- Jackson relation, iFaber fc Jackson! 119761) a rela- 
tionship between stellar m ass and velocity dispersion. 
We adopt the SSP models ofMaraston' ('2005') for Kroupa 
IMF (note that a Salpeter IMF yields M/L ratios, hence 
masses, by a factor 1.6 higher). The M/L ratios (and 
masses) derived in this way from the high density and 
low density relationships vary by only 25 per cent, be- 
cause the younger ages of the low density objects are bal- 
anced by higher metallicities. In the following equation 
we therefore provide the resulting relationship between 
stellar mass and logcr averaged over both environ- 
ments. 

log Af, « 0.63 + 4.52 loga (2) 

This yields the following relations between the stellar 
population parameters and galaxy stellar mass. 

[a/Fe] = -0.459 (-0.459) + 0.062 (0.062) log M^/M^?,) 
[Z/H] = -1.137 (-1.109) + 0.1217 (0.1261) log M*/Mq 
log i/Gyr = 0.427 (0.125) + 0.053 (0.071) log Af^/Mg 

4.1.2. Stellar population parameters 

Fig.|n|shows velocity dispersion and stellar mass versus 
the stellar population parameters age, metallicity [Z/H], 
and [a/Fe] ratio derived from the line indices shown in 
Fig. They are summarized in Table IA3I Observa- 
tional data are the colored symbols (see Fig. EJ, the 
Monte Carlo realizations are shown as black circles. The 
lines indicate the linear relationships of the stellar popu- 
lation parameters with log a used in the simulations (see 
Eqn.^. The residuals of the data and Monte Carlo sim- 
ulations from these relationships are shown in Fig.|7|by 
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Fig. 6. — Stellar population parameters as a functions of veloc- 
ity dispersion (measured within 1/10 re) and stellar mass (upper 
X-axis) for the 'old' subpopulation (circles in Fig. HI and Fig. 151. 
Ages and element abundances are derived with the [a/Fe] enhanced 
SSP models described in Section l2l from the indices H/3, Mgb, and 
(Fe) = (Fe5270 + Fe5335)/2. Colored symbols are the values de- 
rived from the observational data (red: elliptical, blue: SO, green: 
cD), median l-tr error bars are shown. The Monte Carlo realiza- 
tions are shown by the black circles. The latter take the obser- 
vational errors (see Fig. |2J into account plus an intrinsic scatter 
A[a/Fe] = 0.05 dex, A[Z/H] = 0.08 dex, and At = 0.2 ■ t Gyr in 
high density environments (left panel) and Afa/Fe] = 0.07 dex, 
A[Z/H] = 0.10 dex, and At = 0.25 ■ t Gyr in low-density environ- 
ments (right panel). Typical error-bars are given in the bottom- 
right corners. 



the grey and hatched histograms, respectively. Note that 
the simulations take both the observational errors plus 
an intrinsic scatter of the stellar population parameters 
at a given a (A [a/Fe] = 0.05 dex, A[Z/H] = 0.08 dex, 
and At — 0.2 • t Gyr in high density environments 
and A [a/Fe] = 0.07 dex, A[Z/H] = 0.10 dex, and 
Ai — 0.25 • t Gyr in low-density environments) into ac- 
count. The red histograms are simulations without in- 
trinsic scatter for comparison. 

a/Fe We find a very clear correlation between 
a/Fe ratio and logcr or stellar mass. The existence 
of this relation has already been anticipated qual- 
itativ ely by a num ber of authors (e.^. Fisher et al.l 

1 r»nrl J 1 nnr?. t 1 nr\rv'. j 1 o^^^^ 



1995'; Greggio 1997: J0rgcnscn 1999; Kuntschncr 2000^ 
Terlevich fc Forbes 2002: Worthev et al...l992,) . and has 



only recently been confirmed quantitatively by means of 
stellar population mod els wi t h var i able element abun- 
dance ratios_ ( Mehlcrt et al^ |2003t IPro ctor & Sansom 
HOO^; "Thomas ct al. 2002; Tragcr et al. 2000a). The ob- 
served scatter about the relation can only be reproduced 
assuming an intrinsic scatter in [a/Fe] ratio of the order 



0.05 dex as shown by the hatched histograms in Fig.|7| In 
the simulations without intrinsic scatter (red histograms 
in Fig.O, the scatter in a/Fe is significantly smaller than 
the observed distribution. 

Most interestingly, as already shown in 'Thoma s et all 
(2002), the a/Fe-CT relation (and its scatter) of early- type 
galaxies is independent of the environmental density (see 
also Eqn. This reinforc es and exte nds the results 
of t he studies bvlJ0rgensen. F ranx. & Kia ergaardI l|1995D 
and iTerlevich fc Fbrbesr|2002f) . which are based on stel- 
lar population models tied to solar element ratios. Cu- 
riously, early-type galaxies in low density environments 
seem to exhibit stronger CN absorption fea tures than 
in hi gh density environments ISanchcz-Blaz ouez et all 
|2003)- As a significant fraction of C and N comes from 
int ermediate-mass stars in the mass range 5 < M/Mq < 
8 l)Renzini fc Volilll98l() . this points toward about 10^ 
years longer formation timescales in low density environ- 
ments. A delay well above this value can certainly be 
excluded, as the a/Fe ratio is independent of the en- 
vironment. Hence, if true the above result means that 
there must be an amazingly high degree of fine tuning 
in the formation timescales of early-type galaxies as a 
function of environmental density. 

Z/H Besides a/Fe, also total metallicity Z/K cor- 
relates with (7, hence more massive early- type g alaxies 
are also more metal-rich (see also Greggio| (l997l) . This 
result is sensible in terms of chemical evolution, as the 
deeper potential well of massive objects hampers the de- 
velopment of a galactic wind, which leads to a more com- 
plete chemical processing to higher element abun dances 
(e.g.p. rimoto & Yoshii 1987; Edmmids 1990; Mat teuccil 
Il994j) . The observational scatter is best reproduced as- 
suming about 0.1 dex of intrinsic scatter in total metal- 
licity at a given velocity dispersion (hatched histograms 
in Fig.Ej). Different from the a/Fe ratio, however, simu- 
lations without intrinsic scatter still provide a reasonable 
representation of the data (red histograms in Fig. [71). 

Unlike for the a/Fe ratio, we detect a slight impact 
of the environmental density in the sense that galaxies 
in low density environments appear about 0.05-0.1 dex 
more metal-rich than their counterparts in high density 
environments (see also Eqn.^. This implies that galax- 
ies in high-density environments seem to be more efficient 
in ejecting interstellar material in galactic winds, which 
might be caused by galaxy interactions and/or gas strip- 
ping processes. 



Age As noted in prev ious s tudies iKuntschner et alJ 
120011 ITerlevich fc Forbe^l200llivrger et al.ll2000aj) . age 



seems to be the weakest variable among early-type galax- 
ies, the data show no significant trend between age 
and veloc ity dispersion. Howeve r, it has also been 
shown by iKuntschner et alJ l)200H) that correlated er- 
rors of age and metallicity tend to dilute a correlation 
between age and metallicity, because a slight overesti- 
mation of Z/H leads to an underestimation of the age 
and vice versa. In Fig. we show that the observational 
data are best reproduced by a relatively flat but sig- 
nificant correlation between age and velocity dispersion 
(see soli d lines and Eqn. HI) in agreement with previous 
studies ("Proctor. Forbes. & Bcaslev 2004'' IProctor et alJ 
12004: Proctor fc Sansom..2002) . The scatter of the obser- 
vationally derived age (and of the Balmer index H/?, see 
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Fig. 7. — Residuals of the relations shown in Fig.|S|for the observationally derived data (grey histograms) and the simulations (hatched 
histograms). The latter take the observational errors (see Fig. |3J plus an intrinsic scatter A[a/Fe] = 0.05 dex, A[Z/H] = 0.08 dex, and 
At = 0.2 -t Gyr in high density environments (top panels) and A[a/Fe] = 0.07 dex, A[Z/H] = 0.10 dex, and At = 0.25 -t Gyr in low-density 
environments (bottom panels) into account. The red histograms show the residuals of simulations without intrinsic scatter of the stellar 
population parameters. 



Fig. [51) is relatively well reproduced by the simulations 
without any intrinsic scatter (red histograms in Fig.0)- 
A scatter of the order of 20 per cent at a given a improves 
the fit as shown by the hatched histograms in Fig. [7| 

To conclude, the mean age of early- type galaxies is 
weakly correlated with galaxy mass with about 20 per 
cent intrinsic scatter. From this anti-hierarchical rela- 
tionship we infer that the stellar populations of more 
massive early- type galaxies have formed earlier, which is 
further supported by their high a/Fe ratios. 

Unlike in the case of the a/Fe ratio and total metal- 
licity, we find a significant difference between the ages of 
galaxies in high- and low-density environments. There 
is a clear offset indicating that the stellar populations 
of early-type galaxies in low density environments are 
on average abo ut 2 Gyr younger in agreemen t with 
previ ou s results llCaldwell et alJ l2()03t .Kuntschner et al.l 
200l iPoggianti et alJ l2001bl: iProctor et all l2004t 
Terlevich fc ForbesI l2002t iTrager et al.l l2000ar . This 
relative age difference caused by environmental ef- 
fects is in very good agreement with the predic- 
tions from semi-analytic m odels of galaxy formation 

4.2. Low-mass galaxies - triangles 

The objects of the low-mass tail of the velocity dis- 
persion distribution (Fig. ^) with stellar masses < 
10^°Mo (Eqn. 0) exhibit all relatively high Balmer in- 
dices, indicating light-averaged ages between 2 and 3 Gyr 
(triangles in Fig. We recover these objects with 2- 
component composite stellar population models based on 
the the 'old population' (circles in Fig.0J the properties 
of which are specified in Eqn.^ Extrapolation to the ap- 
propriate a yields ages between 6 and 8 Gyr for the base 
old population. The average ages of the low-mass objects 
in both environments are then matched by contaminating 
this base old population with a 10 per cent contribution 
(in mass) of a young 1.2 Gyr) population (see filled 



triangles in Fig. 0J. As mentioned earlier, this is cer- 
tainly not a unique solution, as the contribution in mass 
and the look-back time of the star formation event are 
highly degenerate and therefore only poorly constrained. 
Alternatively, also the addition of a major young com- 
ponent contributing 70 per cent in mass with an age of 
2 Gyr reproduces the observed line indices. We can con- 
clude, however, that relatively recent star formation in 
low-mass early-type galaxies is required to explain the 
observational data (see also Trager et al. 2000a). 

The model with the minor young component gets some 
support from the a/Fe ratios. The bottom panel of Fig.0] 
shows that also the Mg- and Fe- indices are well repro- 
duced by assuming that the young component is 0.2 dex 
less a/Fe enhanced than the base old population. This 
clearly favors the option of a young component over the 
possible presence of old blue horizontal branch stars as 
an explanation for the relatively strong Balmer absorp- 
tion observed in low-mass galaxies. The metal indices of 
the objects are relatively weak, we needed to assume the 
young component to be about 0.2 dex more metal-rich 
than the base old population. 

4.3. Young intermediate-mass galaxies - squares 

This subpopulation of objects with velocity dispersions 
between 110 and 230 km/s, hence 7 • 10^ < M^/Mq < 
2 • 10^^ (see Eqn. I^J, represents a minor fraction of the 
whole sample, namely ^ 14 per cent containing both 
elliptical and lenticular types. There is no obvious de- 
pendence on environmental density and galaxy type. We 
consider the slight bias toward more lenticulars in high- 
density environments (6 SOs out of 8 early-types in high 
density vs. 4 SOs out of 9 early-types in low density) 
as statistically not significant given the relatively small 
sample size. 

The striking feature of these objects is their relatively 
strong Balmer line index, yielding light-average ages be- 
tween 2 and 4 Gyr. To reproduce their H/3 indices, we 
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model these galaxies with a 2-component model adding 
a minor young component as described in the previous 
paragraph. Different from low-mass galaxies, however, 
are the metallicities of the young subcomponent. To 
match the relatively strong [MgFc]' index observed, ex- 
traordinarily high metallicities for the young component 
are required. The smaller the contribution (in mass) of 
the subcomponent is, the higher its metallicity needs 
to be. The avoidance of unplausibly high metallicities 
([Z/H] > 0.8 dex) sets therefore a lower limit to the 
weight of the subpopulation. Metallicities below that 
limit are ensured, if the young subpopulation contributes 
at least 20 per cent in mass. For this lower limit, an age 
of 1.4 Gyr and metallicities around [Z/H] ^ 0.8 dex are 
required to reproduce the observational data. An up- 
per limit can not be constrained. The high metallicity 
of the young component points toward a much higher 
effective yield during the formation of the young com- 
ponent in these intermediate-mass objects, compared to 
the low-mass objects discussed above. This means that 
chemical evolution and the build-up of metals is less effi- 
cient in low-mass objects, in agreement with the general 
picture that objects with shall ower potential wells de- 
velop galactic winds easier (e.g. lArimoto fc Yoshiilll987t 
IChiosi fc Carrarol2002HEdmundsll99nHMatteu cci''199lv. 

Also different from the case of the low-mass galaxies, 
we had to assume the young population to have the same 
a/Fe ratio like the base old population, in order to repro- 
duce their Mg- and Fe-indices (bottom panel of Fig. 01 
compare to dotted model lines for t = 2 Gyr). This 
spoils the interpretation of the enhanced Balmer line 
strengths in terms of the presence of an intermediate- 
age stellar population. The separation in age between 
the old and young components amounts to several bil- 
lion years, plenty of time for Type la supernovae to en- 
rich the interstellar medium with iron. Forming out of 
this Fe-enhanced material, it will be impossible for the 
young component to reach significantly super-solar a/Fe 
ratios. Even the self-enrichment of the young component 
from Type II supernovae during the burst is not sufficient 
to raise the average a/Fe ratio of the newly born stellar 
p opulations indep endent of the burst timescale as shown 
in lThomas et al.l l(T9991 

The interpretation of the observed strong Balmer line 
indices in terms of recent star formation could be re- 
covered under the rather contrived assumption that the 
Type la supernova ejecta originating from stars of the 
base old population have been lost selectively in galactic 
winds before the formation of the young subpopulation. 
Alternatively, the latter must have formed with a stel- 
lar initial mass function (IMF) significantly flatter than 
Salpeter. A flat IMF slope would naturally reproduce 
the very high metallicity of the young component de- 
rived here. On the other hand, compelling evidences for 
systematic variations o f the IMF slope in starbursts have 
not been found so far ijLeithereill 1 998j) . 

But we have to be cautious here. Relatively high a/Fe 
ratios are needed for the young subpopulation, because in 
the models Mg6 decreases significantly with decreasing 
age at highest metallicity, and it is possible that the fit- 
ting functions are in error in this regime owing to bad in- 
put stellar data. On top of this, the extraordinarily high 
metallicities of the young subcomponent ([2'/H] ~ 0.8) 
require extrapolation well above the highest metallicity 



model ([Z/H] — 0.67), which makes the model indices 
very unreliable. The relatively low Mg6 index of the 
model, which causes the need for the relatively high a/Fe 
ratio, might be an artifact of bad calibration. To con- 
clude, lower values for the young subpopulation cannot 
be excluded, which rehabilitates the possibility of the 
presence of an intermediate-age population. 

Alternatively, the strong Balmer lines might be the 
result of blue horizontal branch stars rather than a young 
population, which we shall discuss in detail below. 

4.4. Young massive galaxies - diamonds 

Finally, the diamonds indicate a subclass of very mas- 
sive early- type galaxies (M* > 10^^ Mq) that have 
strong H/3 and strong [MgFe]' absorption indices, imply- 
ing light-averaged ages between 2 and 5 Gyr, and very 
high metallicities of [Z/R] > 0.67 dex. From Fig. H it 
can be seen that the occurrence of these objects is re- 
stricted to high-density environments. Moreover, they 
are only lenticular galaxies, no ellipticals with such low 
average ages and high metallicities are found (see Fig. Q . 
The vast majority (80 per cent) of all SO galaxies with 
(7 > 200 km/s in our sample and the two Coma cD galax- 
ies belong to this object class (see al so Kuntschner 2003 
IMehlert et a,1 ."200,?: Wgia nti et "^12001 bt) . 

The metallicities are more than 0.2 dex above the av- 
erage metallicity of the most massive objects in the old 
base population constrained by Eqn. ^ Such a signifi- 
cant increase in average metallicity can obviously not be 
obtained in a 2-component model with only 10-20 per 
cent of a metal-rich subpopulation, unless unreasonable 
high metallicities above ten times solar are assumed. We 
therefore have to increase the weight of the young pop- 
ulation. The filled diamonds shown in Fig. ^ are models 
in which the old base population is 'perturbed' with a 
2.5 Gyr young population with metallicities as high as 
[Z/H] ~ 0.7 dex contributing 70 per cent in mass. Lower 
metallicities can be adopted only under the assumption 
that an even larger fraction of the stellar populations 
belongs to the second component. 

Interestingly, like in the previous case these objects 
have very high a/Fe ratios. For the young component we 
need to adopt a/Fe ratios as they are given by Eqn. ^ 
hence [a/Fe] > 0.2 for a > 200 km/s. This fits well 
into the picture that the young average ages are not the 
result of minor recent star formation, but rather the con- 
sequence of a recent very massive star formation event 
at a redshift around z = 0.3 forming the majority of the 
galaxy's stellar populations on a short timescale. This 
would imply that massive SOs and cD galaxies in clus- 
ters form very recently, while ellipticals of comparable 
mass contain only old stellar populations. This expla- 
nation seems contrived, and it remains to be seen if it 
can be understood in the framework of galaxy formation 
models. 

4.5. Blue horizontal branch stars at high metallicities 

Alternatively, the strong H/3 indices discussed above 
might be caused by the presence of metal-rich, blue hor- 
izontal branch stars (see Introduction). Such blue hor- 
izontal branch stars may originate from either a metal- 
poor s ubpopulation as discussed in Maraston fc ThomasI 
l)2000() or metal-rich populations with blue horizontal 
branch morphologies due to enhanced mass loss along 
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the red giant branch evolutionary phase. Here we con- 
sider this latter option. 

4.5.1. The origin of blue horizontal branches 

The mass-loss along the red giant branch can be theo- 
retically understood as stellar winds being produced by 
the acoustic flux generated by the superadiabatic con- 
vection in the out er convcctive zone of the red giant 
l)Fusi-Pecci fc R.enzini 1975') . The fraction of acoustic en- 
ergy which is used to remove the stellar envelope is not 
known from physically motivated, first principles. An al- 
ternative is the empirical deter mination of the mass loss 
rate, which has been done by 'Reimers (1975") studying 
the behavior of the circumstellar calcium K line in a large 
number of red giants and supergiants. The normalization 
of the resulting formula, however, suffers from very large 
u ncertainties. 

'Fu si-Pecci fc Renzinil l)1976j) combined the two ap- 
proaches and introduced an efficiency factor rj for 
Reimer's formula. They calibrated the mass loss rate 
with observational constraints like the horizontal branch 
morphology (Rood 197^, the limiting initial mass for 
carbon ignition ('-^ 4 Mq for the Hvades fvan den Heuvell 
[1975 ') ■ and finally the solar mass loss rate. An effi- 
ciency 77 « 0.33 turned out to be the most appropri- 
ate choice. This val ue is adopted in the models of 
IMarastonI l|1998l 120051) . Note that the observational un- 
certainties in Reimer's formula allow this efhciency fac- 
tor to lie in the range 0.3 ^ ?7 ^ 3. The problem is 
that the calibration has been carried out at sub-solar 
metallicity, and deviations from this canonical value, in 
particular at highe r metallicities, cannot be excluded 
i|Greggio fc RenzinillT990|) . 

4.5.2. Effect on the line indices 

The models of TMB are based on IMarastonI ljl998l 
1^05), thus rely on the calibrated value rj = 0.33. Here 
we additionally compute models with the slightly larger 
efficiency factor rj — 0.55, which is well within the range 
of uncertainty allowed by the Reimer's formula. These 
models have bluer horizontal branch morphologies, and 
therefore predict significantly stronger Balmer indices, 
and slightly weaker metallic indices. This is shown in 
the left panel of Fig. |Slby the arrow labeled 'HB' for a 
10 Gyr old population. The arrow labeled 'age' shows 
the location of a model with age t = 1.5 Gyr, which 
yields the same H/3 value, for comparison. The steeper 
slope of the 'HB' model illustrates the weaker sensitivity 
of the metallic line indices to HB morphology than to 
age. The relative effect of Mg b and (Fe) is shown in the 
right panel of Fig. |S1 Relative to (Fe) the Mg b index 
responds stronger to a decrease of age than to a blueing 
of the horizontal branch. As a consequence, the slope of 
the 'HB' model is steeper also in the Mg&-(Fe) diagram. 
Note that the slope is comparable to the model varying 
metallicity at a constant a/Fe ratio and age. 

4.5.3. 2-component model 

These results will mainly affect the age determination. 
To test this, we have constructed a 2-component model 
in which the base old population is contaminated with 50 
per cent of an old metal-rich population with blue hori- 
zontal branch stars. The location of both the squares and 



the diamonds in Fig.0]is reached by these models, hence 
without invoking the presence of intermediate-age popu- 
lations. We had to assume the metallicity of the second 
component to be about 0.2-0.4 dex higher than the av- 
erage of the base population, while the a/Fe ratio needs 
to be about 0.1 dex lower. Hence, the blue HB option re- 
quires the assumption that preferentially the metal-rich 
stars in a composite stellar population develop blue hor- 
izontal branches. The slightly lower a/Fe ratios would 
fit into this picture, as the most metal-rich populations 
must have formed latest and should therefore be least 
a/Fe enhanced. 

4.5.4. Young vs. blue horizontal branch 

To conclude, none of the two options — intermediate- 
age population or blue horizontal branch — can be ex- 
cluded, further observational constraints to distinguish 
the two alternative explanations are required. 

If intermediate-age populations are present in these 
galaxies, they should be detectable by their very red 
near-infrared colors, because of the major contribution 
to the bolometric from asy mptotic giant branch stars at 
ages between 0.5 and 2 Gyr (|Ma,rastonll1 998l I200I I2OO5I: 
IMaTa,stnTi et a Tl l2001l) . 

Un-ambiguous indicators of blue horizontal branches, 
ins tead, s till ne ed to be found. The Call index defined 
bv lR,osel l)1984D . even though promising, does not solve 
the problem as it is still degenerate betwe en gravity and 
metallicity effects. Schiavo n et al.l (|200^ show that in 
metal-poor globular clusters the higher-order Balmer line 
index HSp is more sensitive to the presence of blue hori- 
zontal branches than H/3. Hence, the index ratio HJp /H/3 
may be a tool to break the degeneracy between age and 
horizontal branch morphology. It is not clear, however, if 
this effect is still present at solar metallicities and above, 
which are relevant t o massive galaxies. Moreover, in 
IThomas et all l|200l we show that at such high metal- 
licities the HJp /H/3 index ratio becomes very sensitive to 
the a/Fe ratio in the sense that HSp/HP increases with 
a/Fe. If blue horizontal branches ought to be detected 
in early-type galaxies by means of the H(5f/H/3 ratio, the 
dependence on a/Fe certainly has to be taken into ac- 
count. It will be interesting to make these checks for the 
galaxy sample discussed here. Such an analysis will be 
subject of a future study. 

5. THE ORIGIN OF THE MG-a RELATION 

The correlations of the stellar population parameters 
age, total metallicity, and a/Fe ratio with central ve- 
locity dispersion shown in Fig.|Bland Eqn. fallow us to 
constrain the origin of the well-known Mg-cr relation (top 
panels of Fig. EJ. We recall that the correlations found 
here suggest that all three parameters correlate with cr, 
hence support the Mg-cr relationship. 

Fig. ini shows the observed Mgb-a relation of the whole 
sample high density and low density environments com- 
bined. It can be verified from Fig. |5l that the envi- 
ronment has a negligible effect on the Mg-cr relation 
in agreement with previous findirigs based on much 
larger samples fBernardi et al. | ll998HColless et aDll999t 
Worthey & CoUobcrt 2003). Circles are the 'old pop- 
ulation' representing the major fraction of the sample, 
triangles, squares, and diamonds are the 'young outliers' 
discussed in the previous section (see also Fig. 0J. Red 



12 D. THOMAS, C. MARASTON, R. BENDER, C. MENDES DE OLIVEIRA 




Fig. 8. — Comparison of the effects of age and blue horizontal branch on the absorption indices H/3,Mgf), and (Fe). Symbols and model 
grid like in Fig.|l] The arrows labeled 'HB' show the index changes of a 10 Gyr population ([Z/H] = 0.35, [a/Fe] = 0.3, RGB mass loss 
parameter t] = 0.33) when a blue horizontal branch rj = 0.55 is considered. For comparison, the effect of reducing the age to 1.5 Gyr is 
indicated by the arrows labeled 'age'. 
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Fig. 9. — Mg6-(T relation of the whole sample high density and 
low density environments combined. Colored symbols are observa- 
tional data (red: elliptical, blue: SO, green: cD), median l-cr error 
bars are shown. Black symbols Monte Carlo simulations (symbol 
types like in Fig.|l}. The solid line is the Mg fe-cr relation resulting 
from the relationships given in Eqn. ^averaging over high and low 
density environments. The short-dashed line shows the resulting 
Mg b-a relation if only age is assumed to be the driving parameter 
adopting Eqn. The long-dashed line demonstrates the effect of 
additionally adding the dependence on «/Fe ratio from Eqn. [Tl 



and blue symbols are elliptical and lenticular galaxies. 
Green symbols are cD galaxies. The black symbols are 
the Monte Carlo realizations. The solid line is the Mg b- 
a relation that results from the relationships given in 
Eqn. n (averaging over both environments). The short- 
dashed line shows the resulting Mg 6-cr relation if only age 
is assumed to be the driving parameter adopting Eqn. ^ 
The dependencies of both total metallicity and a/Fe ra- 
tio with a are suppressed. The long-dashed line demon- 
strates the effect of additionally adding the dependence 
on a/Fe ratio from Eqn. 

It can be seen that total metallicity has by far the 
largest share (60 per cent) in shaping the relation, fol- 
lowed by the a/Fe ratio contributing only 23 per cent to 
the observed slope. The contribution from age is about 
17 per cent, hence is smallest but comparable to the one 



of a/Fe. We emphasize that these numbers do not vary 
significantly with environmental density. As can be in- 
ferred from Figs. El and El the scatter of the Mg h-a re- 
lation among the old population (circles) results from 
intrinsic scatter in all three stellar population parame- 
ters. The (Fe)-cr relation is much shallower (see Fig. EJ 
because the increasing a/Fe ratio with a diminishes the 
{Fe)-index, compensating lar gely the increase du e to in- 
creasing metalhcity (see also lMehlert et al.ll2003l) . 

This result fits very well to the fact that also the 
color-magnitude relation appears to be driven by 
meta lli city rather t h an a ge llKodama fc Arimotd 
T997t IMerluzzi et ahl IStanford et al.1 ITflM 

Terlevich et al.l Il999t iVazdekis et alJ I2001D . even 
though the effect of the a/Fe rati o variations on colors 
is sti ll not very well understood ijThomas fc MarastonI 
I2OOI. 

Finally, we briefly examine the location of the 'young 
outliers' of Fig. ^ (triangles, squares and diamonds). 
The low-mass and intermediate-mass objects with young 
light-averaged ages (triangles and squares) are clearly off 
the relation and exhibit systematically lower index val- 
ues (by ~ 0.5 A) at a given o. Their slightly larger light- 
averaged metallicities (see Fig. do obviously not suf- 
fice to bring them back on the relationship. Different the 
class of massive, apparently young Coma cluster lentic- 
ulars (diamonds). Their extremely high light-averaged 
metallicities compensate totally for the decrease in Mg h 
caused by the young light-averaged ages. As far as the 
Mg5-(T relation is concerned, these objects are indistin- 
guishable from the main body of the sample. 

6. FORMATION EPOCHS 

The relations shown in Fig. (HI (see also Eqns.n]and|3Il 

can be used to constrain the epochs of the main star 
formation episodes and the star formation timescales for 
early-type galaxies as a function of their mass and en- 
vironmental densities. It should be emphasized that the 
majority of the objects of the present sample obeys these 
relationships. Only 15 per cent of the sample (squares in 
Fig.^J requires a 20 per cent contribution from interme- 
diate age populations on top of the base old population 
described by Eqn. ^ and further 9 per cent (diamonds) 
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Fig. 10. — Star formation histories of early-type galaxies as a function of their stellar masses M« (see labels). The stellar masses M* 
correspond to the velocity dispersions a = 100, 140, 190, 240, 280, 320 km/s (see Egnl^. The star formation histories are derived from 
the mean ages and [o/Fe] ratios shown in Fig. El and Eqns. □(org) and|ll Redshifts assume = 0.3, = 0.7, Ho = 75 km/s/Mpc 
JSpergel et al.l2003l) . The dotted lines marks the average age of a high density objects with M, = 10^^ Mq {a = 320 km/s) for comparison. 



are best modeled assuming 70 per cent young popula- 
tions. The young low- mass galaxies (triangles in Fig. 0J 
do not play a significant role in the total mass budget 
of early-type galaxies, the latter being dominated by L* 
galaxies with M^, ~ 10^^ Mq situa ted at the turnover 
of the luminosi ty (mass) function ijBaldrv et al.1 1200l: 
IBell et al.ll2003|) . This imphes that about 90 per cent 
of the total stellar mass hosted by early-type galaxies is 
well represented by the relationships given in Eqn. ^ 

6.1. Relating timescale and a/Fe ratio 

The key point of this analysis is that the a/Fe ratio 
lifts the degeneracy between formation epoch and for- 
mation timescale. This is done by relating the observa- 
tionally derived luminosity-averaged a/Fe ratios of the 
stellar populations with timescales of their formation. 

For this purpose we simulate the formation and chemi- 
cal evolution of a set of toy galaxies assuming various star 
formation histories characterized by different timescale s 
similar to the exercise presented in lThomas eraLlljlOOgj) . 

We start with a gas cloud of primordial chemical com- 
position. By imposing a star formation rate as a func- 
tion of time, which is given by the specific star forma- 
tion history assumed, we let the gas cloud form stars. 
These stars enrich the gas of the cloud (the interstellar 
medium) with heavy elements produced mainly in super- 
nova explosions. For the aim of this paper we follow the 
chemical enrichment of a and Fe-peak elements. A key 
point is that we take the delayed enrichment of Fe from 



T ype la supernovae into acc oun t using the prescript ion 
of iGreggio fc Renzinil l)1983l see lThomas et al.lll998l for 

more details). Stellar populations successively form out 
of more and more chemically enriched gas, resulting in 
the build-up of a composite stellar population. Every 
generation of stars in this composite will have the chemi- 
cal composition, hence a/Fe ratio, of the gas out of which 
it had formed. We examine this population 12 Gyr after 
the formation of the first stellar generation, and calculate 
its light averaged [a/Fe] ratio, which corresponds to 
the quantity observationally derived in this study. The 
weight in 1^-luminosity of the single stellar populations 
as a f unction of their a ges is adopted from the M/L ra- 
tios of lMarasto n (2004). The IMF slope is kept fixed to 
Salpeter. Like in Thomas et aL ( 1999 .) , the simulations 
are carried out assuming a closed box, which implies that 
selective mass-loss is assumed not to occur and thus not 
to affect the a/Fe ratio. 

The star formation histories in the different runs are 
chosen to be Gaussians of various widths At (FWHM). 
We find the following linear relationship between a/Fe 
ratio and log At: 

[a/Fe] « 1/5 - 1/6 log At . (4) 

The equation reflects what has already been shown in 
Thomas ct al. (1999), namely that a ratio [a/Fe] = 0.2 
of a composite stellar population requires formation 
timescales At < 1 Gyr. The larger At, the lower is 
the final a/Fe ratio of the stellar population because of 
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the late enrichment of Fe from Type la supernovae. Star 
formation extended over the Hubble time yields an a /Ye 
close to solar. We note that this link between a /Ye ratio 
and star formation timescale is in very good agreement 
with the more d etailed chemical evolution simulations of 
iPipino & Matteucci (^04). 

The combination of Eqns. ^ and 13 with Eqn. ^ then 
yields correlations of At with a and mass for early-type 
galaxies valid in both high and low environmental densi- 
ties. 

logAt«3.44 - 1.68 logo- (5) 

« 3.67 -0.37 logAf*/MQ 

Eqns. n and 13 further give the light-average ages as a 
function of ct, constraining the epoch at which the object 
has formed. 

6.2. The star formation histories 

The resulting star formation histories for various 
galaxy masses and velocity dispersions are shown in 
Fig. E]l The x-axis is lookback time, the upper x- 
axis gives the corresponding redshifts ass uming VLm. = 
0.3, p .\ = 0.7, and Hq = 75 km/s/Mpc ijSDergel et all 
|2003|) . The top and bottom panels show the results for 
dense and loose environments, respectively. 

It should be kept in mind that Fig. ^|is based on av- 
erage ages and a /Ye ratios at a given mass. As discussed 
in Sectional the data require a scatter in a /Ye ratio 
of the order 0.05 — 0.07 dex (see Fig. Cj), which trans- 
lates into a spread in star formation timescale of about 
a factor 2. The observed spread in average age is con- 
sistent with the a scatter of about 20-25 per cent (see 
Fig- Ell- For a cluster early- type galaxy with stellar mass 
A'h ~ 5 • 10" A/q (ct = 280 km/s), for instance, this 
implies a spread about the typical formation redshift at 
2^3 between z ~ 2 and z ^ 5. The blurring caused by 
these effects should be kept in mind for the interpreta- 
tion of the diagrams in Fig.^| which are meant to be a 
schematic illustrations of the mean formation epochs. 

6.3. Dependence on galaxy mass 

The plots summarize what has been already dis- 
cussed in the previous sections. The more massive 
the galaxy is, the f aster and earlier form its stellar 
populations. Indeed iPipino fc Matteuccil l|2004D show 
that the observed relation between a /Ye and a can 
be best reproduced by a model in which gas accretion 
and star formation timescales decrease with increasing 
galaxy mass. Only low-mass galaxies exhibit signifi- 
cant star formation at redshifts below z < 1 as also 
concluded by Brinchmann & Ellis ( 2000) and more re- 
cently ||D£jjii£i2^.£La]J l)2004t) . These results are also in 
very goo d agreement w i th the star formation histories 
found by iGavazzi et al.l l)2002() based on spectrophoto- 
metric data. Interestingly, the link between star forma- 
tion histories and galaxy mass shown in Fig. 1101 agrees 
well wit h the predictions from N- body-tree-SPH simula- 
tions bv lChiosi Ca,rra,rnl lf20(B . which are based on a 
monolithic collapse inside a dark matter halo rather than 
hierarchical structure formation. 

Note that we have assumed continuous star formation 
in order to keep the number of free parameters low. The 
principal aim of this exercise is to illustrate the typi- 
cal timescale of the star formation histories in galaxies 



as a function of their masses. In reality, star forma- 
tion in individual objects will be less smooth and will 
be shaped by several bu rsty star formation episodes (see 
iChiosi fc Garrard I2002|) . The functions of Fig. [TUl thus 
have to be understood as integrals at a given velocity dis- 
persion over bursty star formation histories. In particular 
the extended star formation histories of the low-mass ob- 
jects might rather be understood in terms of episodically 
returning star forming bursts, simila r to what is found 
in dwarf galaxies of the Local Group (Mateo 1998). 

In massive systems, the burst might well be triggered 
by merger events. Indeed, a considerable fraction of 
massive ellipticals have kincmatically decoupled cores 
|BendeiE988 ; Franx & lUingworth 1988), the presence 
of which requires the merger to be accompanied by dissi- 
pational processes and star formation. Interestingly, for a 
number of ellipticals with kincmatically decoupled cores, 
it has been found that their stellar populations are ho- 
mogeneously o ld an d g/Fe enhanced (iDavies et al.ll2001t 
iMehlert et a ll 119981 iSurma fc Bended 1199,'iD . indicating 
the merger to have happened at early times, in agree- 
ment with the formati on histories of Fig. [TUl (a counter- 
example is NGC 2865. iHau. Garter, fc Balcell!3[T99l . 

6.4. Influence of the environment 

The second main result is that galaxy formation 
appears delayed by ~ 2 Gyr in low d ensity environ- 
ment s c onfirming previou s findin gs iC^^we^^^^ 
20031: iKuntschner et all 120021: iPoggianti et al 
aOQlbt iProctor et alJ 120041: iTerlevich fc YorheA 12002 
■Trager et al. 2000a). 

This is in agreement with studies of the field 
galaxy pop ulation at intermediat e redshifts a round 
z ^ 1 llGebhardt et all l2003t iTreu et alJ l200l 
Ivan Dokkum et alJl200lD . While the formation of mas- 
sive early-type galaxies in clusters sets in at redshifts 
around z ^ 5 and is completed (in terms of their stel- 
lar populations) at about z ~ 2, a considerable fraction 
of massive galaxy formation in low density environments 
occurs at lower redshifts between z ~ 1 and z ~ 2. We 
recall here that the a/Ye-a relation does not depend on 
environmental density. The timescale of star formation 
in an object is determined by its (final) mass rather than 
the environment. As a consequence, massive early-type 
galaxies in low density environments form later but on 
similar timescales like their counterparts in high environ- 
mental densities as illustrated in Fig. EH 

7. DISCUSSION 

7.1. The local abundance pattern at high metallicities 

In a recent study, IProctor et al.l ((203) argue that the 
correlation of a/Fe with velocity dispersion is an 'arte- 
fact of incomplete calibration (of the models) to the Lick 
system '. The basis fo r their conclusion is the recent find- 
ing of iBensbv eTall (Hooi), that the Milky Way disk 
stars with super-solar metallicities do not have solar but 
sub-solar 0/F e ratios, 0/Fe decreasing with increasing 
Fe/H (see also lAUende Prieto et al.ll2004HPompeia et all 
I2003D . Stellar population models, being calibrated with 
these stars, are then likely to reflect sub-solar 0/Fe ra- 
tios at high metallicities. If the element Oxygen is taken 
as representative for all a-elements, this would obviously 
lead to an overestimation of the 0/Fe ratio (or a /Ye ra- 
tio) in stellar populations with super-solar metallicities 
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like early-type galaxies. iProctor et al.l l)2004|) show that 
the correction for this effect makes the a/Fe-cr relation 
disappear. 

However, there is a further v ery important (and s ome- 
what puzzling) point shown in' Bensbv et alJ l)2004f ^. All 
the other a-elements, Mg, Si, Al, Na, Ca, and Ti, do 
not show this pattern, but level-off at a s olar X/Fe ra- 
tio in stars with super-solar metallicities ijBensbv et al.l 
|2003; Fuhrmann 1998) . Only the element oxygen (and 
ma ybe also carbon, iThoren. Edvardsson. fc Gustafssonl 
l2004i|') sticks out of this homogeneous behavior. This 
differential pattern between Oxygen and the other a- 
elements cannot be neglected in models taking the lo- 
cal abundance pattern at high metallicities into account. 
This is particularly important in the context of this work, 
as the a/Fe ratios presented here are derived from the 
absorption line indices Mg b and (Fe) , which are sensitive 
mostly to the element s Fe an d Mg, but not O ijKorn et al.l 
EdOS; Tripicco & Rdi rmol . 

We investigated this issue constructing two additional 
flavors of stellar population models that account for this 
local abundance pattern at high metallicities. The first 
model takes oxygen as rep resentative for all a-elements 
(like in IProctor et al.ll2004|) . in the second model oxygen 
is detached from the rest of the a-elements as discussed 
above. The resulting models in the index-index planes 
and the stellar population param eter s deri ved with these 
models are show n in F igs. IBIII and IB 131 (equivalent to 
Fig. El and Figs. IbT2I and IbTI feauivalent to Fig.EJ in 
the Appendix. The results are discussed in detail there. 

To summa r ize, we reproduce the result of 
IProctor et al.l l)2004f) with the model in which oxy- 
gen and a-element are not separated. The a/Fe-a 
relation disappears, the Z/Fl-a relation becomes tighter, 
while the ages remain unchanged. When oxygen is 
detached from the other a-elements as suggested by ob- 
servations, the picture is entirely different as explained 
in detail in Section ^ The correction for the local 
abundance pattern has no significant impact on the 
model, when the non-solar 0/Mg ratios of metal-rich 
stars in the solar neighborhood are taken into account. 
The model basically moves back to its original position. 
As a consequence, the stellar population parameters 
derived with this model are inexcellent agreement with 
the ones presented in this work (see Fig. lST4|l . 

To conclude, the results of this paper are robust against 
the correction for the local abundance pattern at high 
metallicities. In particular, the a/Fe-a relation found 
here is real and is not caused by calibration artefacts. 
What remains to be understood is the discrepancy be- 
tween oxygen and the other a-elements in the metal-rich 
stars of the solar n eighborhood. 

As discussed in iBensbv et all ll2004|) . a possible ex- 
planation could be that Type la supernovae do pro- 
duce a significant amount of magnesium, which keeps 
the Mg/Fe ratio solar, while 0/Fe decreases with in- 
creasing metallicity. This would also fit to the slightly 
higher 0/Fe than Mg/Fe ratios i n metal-poor stars of the 
Milky Way ijBensbv et alJl200-f) . In this case, the stan- 
dard model for Type la nucleosynthesis (the so -called 
'W7' model; iNomoto. Thielemann, fc YokoillTosl needs 
to be revised. This leads to the speculation that the 
importance of Type la enrichment may increase steadily 



with atomic number of the element up to Fe and Ni. If 
this is the case, the higher the atomic number of the el- 
ement, the flatter should be the X/Fe-Fe/H relation in 
the Milky Way. In particular elements with high atomic 
number should have lower X/Fe ratios in the metal-poor 
halo stars, wh i ch ha s not been confirmed in the review 
of McWillia^ l(T99l . Maybe new data sets will revo- 
lutionize this statement. If true, the low Ca/Mg ratios 
of early-type galaxies (e.g. JThomas et al. 2003a) would 
naturally be explain ed in terms of de layed Type la su- 
pernova enrichment ijPipino fc Matteucci 200 4'). 

The relevant inference for this paper would be that the 
formation timescales discussed in Section El (Eqn. Q are 
underestimated. A detailed investigation and quantifica- 
tion of this effect goes far beyond the scope of this paper, 
and in particular cannot be performed in a straightfor- 
ward way as the appropriate Type la nucleosynthesis cal- 
culations are not available. Note however, that the high 
a/Fe ratios derived here for massive galaxies, suggesting 
short star formation timescales, are coupled with old av- 
erage ages, which further reinforces this interpretation. 
The results of this work suggest star formation timescales 
as short as only a few 10® years in the most massive early- 
type galaxies with M* ~ IQ-*^^ Mq. It would certainly be 
desirable if such extreme timescales were slightly relaxed 
by the additional contribution from Type la supernovae 
to the enrichment of the heavier a-elements. It should be 
emphasized that a fundamental conclusion of this paper 
is that low-mass galaxies have more extended formation 
timescales. This purely differential result is not affected, 
if Type la supernovae represent an additional source of 
a-elements. 

7.2. Comparison to early results 

The findings discussed in this paper are supported by 
a number of evidences from the tightness and the red- 
shift evol ution of the scaling relations of ea r ly-type galax- 
ies (e.g. irrag6n-Sa,]a,ma,ricaet all 179931 IRender et alJ 
[l9M IBower et all 119921: IKodama et all 119991: 
'Renzini fc Ciotti"1993': 'Sag ha et alJ l2000': 'Ziegler et alJ 
2001; Ziegler ct al. 1999, see review by Peebles 2002). 
It is known since these early (and some more recent) 
works that the formation ages of the stellar populations 
in elliptical galaxies have to be above redshifts z > 2. 
New in this paper is the detailed disentanglement of star 
formation histories as a function of galaxy mass and the 
distinction between dense and loose environments. 

We confirm such high formation redshifts for 
intermediate-mass and massive (Af* > 10^*^ Mq 
or tJ > 140 km/s) early- type galaxies in dense 
environments in good agreement also with the re- 
cently discovered relationshi p between globular cluster 
color and galaxy luminosity l)Strader. Brodie. fc Forbed 
120041) and the K-z relation o f powerful radio galaxies 
l)R,occa-Volmerange et al.l|200 4>) . 

We find, however, that ellipticals in looser environ- 
ments form at somewhat lower redshifts between z ~ 2 
and z 1. Low-mass galaxies (M* < 2 ■ 10^" Mq or 
a < 140 km/s) in low density regions apparently suffer 
from significant star formation episodes below redshift 
z ^ 1. This paper further adds the conclusion that the 
stellar populations of more massive early-type galaxies 
form earlier and faster. This result con f irms e arly, more 
qualitative estimates by iBender et alJ l)1996(l based on 
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the redshift evolution of the Mg-a relation, who had al- 
ready noted that the formation redshifts of the most mas- 
sive early-type galaxies in clusters may be even above 
z > 4. Here we find that while the formation epoch 
is delayed in low-density environments by about 2 Gyr, 
the timescale is independent of the environmental density 
and depends only on the galaxy's mass concentration. 

7.3. The stellar mass density at z > 1 

The above evidences provide consistent constraints on 
the epochs when the stellar populations of early-type 
galaxies have formed. Our results imply that at red- 
shift around z ^ 1, the formation of the total stellar 
mass in early-type galaxies with Af* > 2 • 10^° Mq or 
a > 140 km/s in both high and low density environments 
is almost completed. Stellar mass growth in early-type 
galaxies below redshift z ^ 1 is restricted to low-mass 
objects in low density environments. As the latter make 
up only a small fraction of the total mass of a galaxy 
population, despite them being numerous, we can infer 
that stellar mass growth in early-type galaxies is largely 
completed around redshift z ~ 1. 

This galaxy type is expected to host betwe en 1/2 and 
3/4 of the total ste llar mass in the universe iBell et al.l 
l2nn3HRen zini"199;j^. As also massive disks arc already 
present at redshifts between 1 and 2 ( Ravindranath ct al. 
l2004HStockton et al.ll20'o^ . this imphes that at least 50 
per cent of the total stellar mass density must have al- 
ready been formed by z ^ 1. Our results further sug- 
gest that a significant increase of the stellar mass den- 
sity between redshifts 1 and 2 should be present caused 
mainly by the field galaxy population, while at redshifts 
around 3 the stellar mass density should be dominated 
by galaxies in high density environments. Depending on 
the exact partition between objects in high and low en- 
vironmental densities, this implies the stellar mass den- 
sity around redshift 3 to be as low as 10-30 per cent 
of the local value. Observational estimates of the to- 
tal stellar mass density as a function of redshift are in 
overall good agreement with this result. Recent stud- 
ies are consistent with only mild evolution of the order 
20 — 30 per cent up to z ~ 1 followed by a significant 
drop in ste llar mass d ensity bet ween reds hifts 1 and 2 
(IBell et al.l l2003: Brinc hmann fc Ellis 200 0l lCohenll200l 
[Dickinson ct al. 2003: Fontana et al. 2004: Fontan a et aP 
11)03; Glazcbrook et al. 2004; Rudnick ct al. 2003). Our 
result suggests that the significant mass growth in the 
early-type galaxy population below z ~ 1 observed by 
IBell et all l)2004|) is most likely caused by less massive 
objects. 

7.4. Climbing up the redshift ladder 

One fundamental consequence from the results of this 
paper is the prediction that very massive objects A/» > 
10^^ -^0) with very high star formation rates up to sev- 
eral 1000 Mo/yr should be present between redshifts 2 
and 5. 

Such objects have indeed been recently found. At least 
part of the SCUBA sources at 2 < z < 3 turn out to be 
massive star forming galaxies with extrem ely high star 
form ation rates up to several 1000 Mc^/y v (jGenzel et al.l 
l2003HSmai l ct al. 2002':' Tecza et al.l2004l) . These objects 
are likely to be the precursors of the most massive ellipti- 
cals (M* w 10^^ -^0) forming in a violent star formation 



episode at high redshift. Ly-break galaxies, instead, ex- 
hibit more mode rate star formation rates of the orde r 
a few 10 Mg/yr ijBaker et alJl2004HPettini et alJl200H) . 
and may therefore be the precursors of less massive early- 
type galaxies in Fig. ITT)| (CM^ ^ 10 " Mq) and/or spi- 
ral bulges. Moreover. Pettini et all hWl) ?iii<l cB58 (a 
lensed galaxy at redshift z = 2.73. iSeitz et al.lll998(l 
to have a significantly a-enhanced interstellar medium^ 
which is in good agreement with the a /Fe- enhancement 
derived in this paper for the stellar population of local 
ellipticals. 

A further direct inference from Fig. EH is that 
massive old, passively evolved galaxies should ex- 
ist at intermediate redshifts between 1 and 2. 
Again, these objects have been found. Medium- 
deep infrared surveys detect a relatively large num- 
ber of so-called extremely red objects (EROs) be- 
tween z ~ 1 and 1.5, whose SEDs b e st fit evolved 
stellar pop u lations llCaputi et a.1.1 l2004t IC imatti e t alJ 
1200^, "2004: 'Daddi et al.' ■2000'; "Moustakas et al.' '2001 
Pozzct ti fc Mannucci 2 000; Saracco ct al. 2003, 200j 
Vais anen fc JohanssonI l2004() . The evolved i^T-selected 
EROs mentioned above s how evidence for strong clus- 
tering l)Daddi et al.l 1200!^ . hence are likely candidates 
for being the precursors of the cluster objects shown in 
the top panel of Fig. EH In a recent paper, Pierini et al] 
12004) show that the colors of some 'evolved' EROs are 
also consistent with a dusty posi-starburst population 
with ages around 0.5 Gyr, which fits to the star forma- 
tion histories shown in Fig. 1101 for the low density envi- 
ronment. It would be interesting to investigate the clus- 
tering properties of this subclass of EROs, to verify that 
they preferably are found in low density environments. 
Also in agreement with our result is the recent finding 
in if-selected galaxy surveys that the most luminous ob- 
jects are best fit by older stellar populations at all red- 
shifts up to z 1, which implies the typical M/ L ratio 
(in K) of more m assive galaxies to be larger than that of 
less m assive ones l|Fontana et al.ll200^ iGlazebrook et alJ 

Finally, deeper redshift surveys are now finding 
evolved objects even at redshifts around 2 and 3 
llGa,nnti et a,lJ l2005t IGhen Marzkel IM iFraux et a] J 
120031 iMcCarthv et al.ll2004^ and massive star forming 
objects around redshifts 2 a nd 5 fd e Mellq et alJ 120041: 
iLehnert fc Bremed 120031 |^tani et alJ l200H) upto z ~ 
6.5 l|Kurk e t al. 2004V again in good agreement with the 
prediction from the properties of local galaxies made in 
this paper. 

To conclude, it seems the two approaches, archaeol- 
ogy of stellar populations in local galaxies and redshift 
surveys, are converging to a consistent picture in which 
massive galaxy formation occurs at fairly high redshifts 
accompanied by violent star formation. 

7.5. Models oj galaxy formation 

This consistent picture of high formation redshifts 
of the stellar populations in massive (early-type) 
galaxies severely challenges the co nventional semi- 
analytic mo dels of galaxy formation jCole et all l2000|: 
H^^to^^^^. 2003; Kauffmann et al. 1993; Menci et alJ 
'2003^ iSome'rville fc Primack 1999), and their modifica- 
tion h as become unavoidable. As reviewed bv <Somervill3 
()2004l) . recent renditions of both semi-analytic models 
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and hydro dynamical simulations of hierarchical galaxy 
formation l|Nagamine et al.l2004l) are consistent with the 
stellar mass function around z = 1 and with the cos - 
mic star formation history (see also lFontana et a'l]l2004|) , 
hence global quantities are well reproduced. However, 
the models have difficulty producing the old average ages 
of massive galaxies in the local universe, and accordingly 
the number density of luminous red galaxies at interme- 
diate redshift as well as the high sta r formation rates 
in massive objects at high redshifts ijSomervi 
iSomerville eT^^OO^ . 

At least part of the reason for this shortfall may be con- 
nected to the way non-baryonic dark matter and bary- 
onic matter are linked. In current models, star formation 
is tightly linked to the assembly history of dark mat- 
ter halos, so that galaxies with longer assembly times 
also form stars on longer timescales. A possible solution 
to reconcile extended assembly times with fast star for- 
mation histories is a scenario, in which ma ssive objects 
form their stars early and assemble later ffiow er et al.l 
[1998; van Dokkum ct al. 1999). This picture gets some 
theoretical support from the fact that the kinematics of 
massive early-type galaxies require the origin of merg- 
ers between bulge-dominated rather than disk galax - 
ies lIKhochfar fc BurkertI 1200^ INaab fc BurkertI l20nl . 
However, in this case, it would be difficult for early-type 
galaxies to establish the correlation between metallicity 
and mass reported here. More fundamental modifica- 
tions of the models are certainly necessary. 

iGranato et al.l l|2001l ^ have presented a promising — 
even though quite heuristic — approach, in which star 
formation is enhanced in massive systems because of 
feed back between star forrnation and AGN activity (see 
also Cattaneo fc Bernardil 120031 and iKawata fc Gibsoiil 
20041). This i dea is now put on physical grounds in 
Granato et alJ ((20^, where the evolution of gas in a 
dark matter halo is described as a function of gravity, 
radiative cooling, heating by feedback from supernovae, 
and — most importantly — heating by feedback from the 
growing active nucleus. 

The assumption of a tight interplay between star for- 
mation and nuclear activity gets indeed strong obser- 
vational support from the fact that besides the a/Fe 
ratio, also black hole mass correlates with velocity 
dispersion and galaxy mass (iFerrarese fc Mcrritt 2000; 
iGebhardt et al.ll2000l[i . T he resulting 'anti-hi erarchical 
baryonic collapse' in the IGranato et alJ 1)20041 model 
leads to higher formation redshifts and shorter formation 
timescales of the stellar populations in massive halos. In 
this way bot h the relatio nship between a/Fe ratio and 
galaxy mass l|R,omano et al.ll200'^ . and the presence of 
old galaxies at z ~ 1 (EROs) and of violently star form- 
i ng objects at z ^ 3 (SC UBA sources) are reproduced. 

IGranato et al.l 1)20041) consider the evolution of mas- 
sive halos in the framework of semi-analytic models only 
above redshifts z > 1. Hence, a fundamental problem 
that remains to be solved is the accretion of cold gas 
onto the central galaxy of a dark matter halo at lower 
redshift leading to late star formation and the galaxy's 
rejuvenation. 

An important step to a solution of this la t ter pr oblem 
might have been recently made by .Binnevl l)2004|) , who 
suggests a scenario in which gas cooling is the key pro- 



cess. Unlike small systems with shallow potential wells, 
massive galaxies cannot eject but retain their hot gas, 
the cooling of which is inhibited by an episodically ac- 
tive galactic nucleus. This 'atmosphere' of hot gas builds 
up a shield around the galaxy preventing cool gas to fall 
in and to form stars. The inclusion of this effect in semi- 
analytic models of hierarchical galaxy formation would 
certainly be the right step toward reconciling our cur- 
rent understanding of galaxy formation with the obser- 
vational constraints set by the stellar populations prop- 
erties of massive galaxies. 

8. CONCLUSIONS 

We study the st ellar population properties of 124 
early-type galaxies l)Beuing et alJ |200 2\ Gonzalez '1993"; 
lMehlert~ al. .2003,} with velocity dispersions between 
50 and 350 km/s in both high and low density envi- 
ronments. From the absorption line indices H/3, Mg&, 
and (Fe) we derive ages, metallicities Z/Fi, and a/Fe el- 
ement ratios using our abundance ratio sensitive stellar 
population models (TMB). The a/Fe ratio is the key to 
determine the formation timescales of the objects, as it 
quantifies the relative importance of the chemical enrich- 
ment from Type II and the delayed Type la supernovae. 
Together with average ages, we use this information to 
set constraints on the epochs, when early-type galaxies 
form the bulk of their stars. To get a handle on the 
effect of correlated errors, we compare the observed sam- 
ple with mock galaxy samples with properties as close 
as possible to the observational data in terms of sample 
size, velocity dispersion distribution, and measurement 
errors. 

For the bulk of the sample we find that all three stellar 
population parameters age, total metallicity, and a/Fe 
ratio correlate with velocity dispersion, hence galaxy 
mass. We show that these results are robust against 
recent revisions of the local abundance pattern at high 
metallicities. To recover the observed scatter we need to 
assume an intrinsic scatter of about 20 per cent in age, 
0.08 dex in [Z/R] and 0.05 dex in [a/Fe]. Both zero- 
point and slope of the a/Fe-a relation are independent 
of the environmental density, while the ages of objects 
in low density environments appear systematically lower 
accompanied by slightly higher metallicities. 

We identify three classes of 'outliers' from these rela- 
tionships. 1) The low-mass objects with A/* < 10^" Mq 
(cr < 130 km/s) in our sample show clear evidence 
for the presence of intermediate-age populations. Their 
line indices are best reproduced with a 2-component 
model in which a with a 1.5-2 Gyr population with low 
a/Fe ratios is added to a base population defined by 
the above relationships. 2) About 20 per cent of the 
intermediate-mass objects with 10^° < M^/M© < 10^^ 
(110 < cr/km/s < 230) must have either a young sub- 
population or a blue horizontal branch. Compared to 
the low-mass objects, the case for the intermediate- age 
population is less convincing because of relatively high 
a/Fe ratios. 3) Most intriguing is a population of mas- 
sive (M* > 10" Mq or (t > 200 km/s) SO galaxies, only 
present in the high density sample, that require a major 
fraction of their stars (~ 70 per cent) to have formed 
relatively recently at lookback times of about 2-3 Gyr, 
corresponding to z 0.2. Note that these objects, even 
though being interesting in themselves, represent only a 
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minor fraction of the total stellar mass budget in early- 
type galaxies. 

The increase of total metallicity, indicating the com- 
pleteness of chemical processing, mainly reflects the in- 
ability of massive galaxies with deeper potential wells to 
develop galactic winds. The increase of age and a/Fe 
ratio with a shows that more massive galaxies form their 
stellar populations both earlier and faster. Our results 
suggest that the well-known Mg-tr relation is mainly 
driven by metallicity (~ 60 per cent contribution) with 
almost equal share from a/Fe ratio (~ 23 per cent) and 
age (~ 17 per cent). The increase of a/Fe with velocity 
dispersion is responsible for the much weaker correlation 
of (Fe) with galaxy mass. 

With a simple chemical evolution model that takes the 
delayed enrichment from Type la supernovae into ac- 
count, we translate the ages and a/Fe ratios of the above 
relationships into star formation histories as a functions 
of galaxy's velocity dispersions and environmental den- 
sities. We show that the higher ages and a/Fe ratios of 
more massive early-type galaxies imply faster formation 
timescales at earlier epochs. The lower the mass of a 
galaxy is, the more it suffers from episodes of late star 
formation extending to redshifts below z ~ 1 for galax- 
ies with < 10^" Mq. Massive early-type galaxies 
(M* > 10^^ Mq) in high density environments form their 
stellar populations between redshifts 2 and 5. Galaxies 
in low density environments appear on average 1-2 Gyr 
younger than their counterparts in high environmental 
densities, but obey the same relationship between a/Fe 
ratio and velocity dispersion. This implies that they form 
on the same timescales, but the epochs of their formation 
are delayed to lower redshifts between z ~ 1 and z ^ 1. 
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The 'archaeology approach' pursued in this paper and 
observations of the high-redshift universe are now con- 
verging to a consistent picture. The results of this work 
imply that at least 50 per cent of the total stellar mass 
density must have already been formed by z ^ 1, which 
is in good agreement with observational estimates of the 
total stellar mass density as a function of redshift. We 
conclude that significant mass growth in the early-type 
galaxy population below z ^ 1 must be restricted to less 
massive objects, and a significant increase of the stellar 
mass density between redshifts 1 and 2 should be present 
caused mainly by the field galaxy population. As shown 
here, the properties of local galaxies predict the presence 
of massive objects with very high star formation rates 
of the order 1000 A/0yr~^ around redshifts 2-5, and the 
existence of old evolved galaxies at redshifts between 1 
and 2, both now confirmed observationally. 

Current models of hierarchical galaxy formation are 
consistent with global properties like the stellar mass 
function, but remain challenged by the ages of local mas- 
sive galaxies , the number densities of luminous red ob- 
jects at intermediate redshift, and the high star forma- 
tion rates in massive objects at high redshifts. The prob- 
lem seems tightly linked to the (uncertain) treatment of 
baryonic matter in dark halos, and likely solutions lie 
in the inclusion of AGN feedback and refinement of gas 
heating/cooling processes. 
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APPENDIX 

DATA SAMPLE 

In Table [X2l we list list velocity dispersions, line indices H/3, Mg6, and (Fe), galax y ty pes, environmental densities, 
literature sources, and classifications of Fig. 01 of the sample investigated. Table IA3I gives the stellar population 
parameters age, metallicity, and a/Fe ratio derived with the TMB models. The l-tr errors are estimated through 
Monte Carlo simulations. We perturb the line indices assuming a Gaussian probability distribution defined by the 1-a 
measurement errors quoted above. The resulting l-cr errors correspond to the widths of the distributions in the stellar 
population parameters from 100 Monte Carlo realization per object. 



THE EPOCHS OF EARLY-TYPE GALAXY FORMATION 

Table A2. Sample parameters and line indices 
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Name 


Sigma 


Error 




Error 


Mg6 


Error 


(Fe) 


Error 


Type 


Environment 


Source 


Class 


NGC 0221 


72. 


,1 


2. 


.6 


2 


.31 


0.05 


2.96 


0.03 


2.75 


0.03 


E 


Low density 


Gonzalez 


triangles 


NGC 0224 


156, 


,1 


3 


.7 


1 


.67 


0.07 


4.85 


0.05 


3.10 


0.04 


SO 


Low density 


Gonzalez 


circles 


NGC 0315 


321, 


,0 


3 


.9 


1 


.74 


0.06 


4.84 


0.05 


2.88 


0.05 


E 


Low density 


Gonzalez 


circles 


NGC 0507 


262, 


,2 


6 


,4 


1 


.73 


0.09 


4.52 


0.11 


2.78 


0.10 


E 


Low density 


Gonzalez 


circles 


NGC 0547 


235, 


,6 


3 


.6 


1 


.58 


0.07 


5.02 


0.05 


2.82 


0.05 


E 


Low density 


Gonzalez 


circles 


NGC 0584 


193, 


,2 


2 


.8 


2 


.08 


0.05 


4.33 


0.04 


2.90 


0.03 


SO 


Low density 


Gonzalez 


squares 


NGC 0636 


160, 


,3 


2 


.9 


1 


.89 


0.04 


4.20 


0.04 


3.03 


0.04 


E 


Low density 


Gonzalez 


circles 


NGC 0720 


238, 


,6 


5, 


,2 


1 


.77 


0.12 


5.17 


0.11 


2.87 


0.09 


E 


Low density 


Gonzalez 


circles 


NGC 0821 


188, 


,7 


2 


.9 


1 


.66 


0.04 


4.53 


0.04 


2.95 


0.04 


E 


Low density 


Gonzalez 


circles 


NGC 1453 


286, 


,5 


3 


.5 


1 


.60 


0.06 


4.95 


0.05 


2.98 


0.05 


E 


Low density 


Gonzalez 


circles 


NGC 1600 


314, 


,8 


4, 


.3 


1 


.55 


0.07 


5.13 


0.06 


3.06 


0.06 


E 


Low density 


Gonzalez 


circles 


NGC 1700 


227, 


,3 


3 


.0 


2 


.11 


0.05 


4.15 


0.04 


3.00 


0.04 


E 


Low density 


Gonzalez 


squares 


NGC 2300 


251, 


.8 


3 


.3 


1 


.68 


0.06 


4.98 


0.05 


2.97 


0.05 


E 


Low density 


Gonzalez 


circles 


NGC 2778 


154. 


,4 


3, 


,2 


1 


.77 


0.08 


4.70 


0.06 


2.85 


0.05 


E 


Low density 


Gonzalez 


circles 


NGC 3377 


107, 


.6 


2, 


.6 


2 


.09 


0.05 


3.99 


0.03 


2.61 


0.03 


E 


Low density 


Gonzalez 


squares 


NGC 3379 


203. 


,2 


2, 


.7 


1 


.62 


0.05 


4.78 


0.03 


2.86 


0.03 


E 


Low density 


Gonzalez 


circles 


NGC 3608 


177. 


,7 


3. 


.0 


1 


.69 


0.06 


4.61 


0.04 


2.94 


0.04 


E 


Low density 


Gonzalez 


circles 


NGC 3818 


173. 


,2 


4, 


.2 


1 


.71 


0.08 


4.88 


0.07 


2.97 


0.06 


E 


Low density 


Gonzalez 


circles 


NGC 4261 


288. 


.3 


3 


.3 


1 


.34 


0.06 


5.11 


0.04 


3.01 


0.04 


E 


High density 


Gonzalez 


circles 


NGC 4278 


232. 


,5 


2 


.8 


1 


.56 


0.05 


4.92 


0.04 


2.68 


0.04 


E 


Low density 


Gonzalez 


circles 


NGC 4374 


282. 


,1 


2 


.8 


1 


.51 


0.04 


4.78 


0.03 


2.82 


0.03 


E 


Higii density 


Gonzalez 


circles 


NGC 4472 


279. 


,2 


4, 


,4 


1 


.62 


0.06 


4.85 


0.06 


2.91 


0.05 


E 


High density 


Gonzalez 


circles 


NGC 4478 


127. 


,7 


4, 


.4 


1 


.84 


0.06 


4.33 


0.06 


2.94 


0.05 


E 


High density 


Gonzalez 


circles 


NGC 4489 


47. 


,2 


3 


.8 


2 


.39 


0.07 


3.21 


0.06 


2.66 


0.05 


E 


Higli density 


Gonzalez 


triangles 


NGC 4552 


251. 


,8 


2 


.8 


1 


.47 


0.05 


5.15 


0.03 


2.99 


0.03 


E 


High density 


Gonzalez 


circles 


NGC 4649 


309. 


,8 


3 


.0 


1 


.40 


0.05 


5.33 


0.04 


3.01 


0.03 


SO 


High density 


Gonzalez 


circles 


NGC 4697 


162. 


,4 


3 


.9 


1 


.75 


0.07 


4.08 


0.05 


2.77 


0.04 


E 


High density 


Gonzalez 


circles 


NGC 5638 


154. 


,2 


2 


.8 


1 


.65 


0.04 


4.64 


0.04 


2.84 


0.04 


E 


Low density 


Gonzalez 


circles 


NGC 5812 


200. 


,3 


3 


.0 


1 


.70 


0.04 


4.81 


0.04 


3.06 


0.04 


E 


Low density 


Gonzalez 


circles 


NGC 5813 


204. 


.8 


3 


,2 


1 


.42 


0.07 


4.65 


0.05 


2.67 


0.05 


E 


Low density 


Gonzalez 


circles 


NGC 5831 


160. 


.5 


2 


.7 


2 


.00 


0.05 


4.38 


0.04 


3.05 


0.03 


E 


Low density 


Gonzalez 


diamonds 


NGC 5846 


223. 


.5 


4, 


.0 


1 


.45 


0.07 


4.93 


0.05 


2.86 


0.04 


SO 


Low density 


Gonzalez 


circles 


NGC 6127 


238. 


.9 


4, 


.4 


1 


.50 


0.05 


4.96 


0.06 


2.85 


0.05 


E 


Low density 


Gonzalez 


circles 


NGC 6702 


173. 


.8 


3. 


.0 


2 


.46 


0.06 


3.80 


0.04 


3.00 


0.04 


E 


Low density 


Gonzalez 


squares 


NGC 6703 


182. 


.8 


2. 


.9 


1 


.88 


0.06 


4.30 


0.04 


2.93 


0.04 


SO 


Low density 


Gonzalez 


circles 


NGC 7052 


273. 


.8 


3. 


.6 


1 


.48 


0.07 


5.02 


0.06 


2.84 


0.05 


E 


Low density 


Gonzalez 


circles 


NGC 7454 


106. 


.5 


2. 


.9 


2, 


.15 


0.06 


3.27 


0.05 


2.48 


0.04 


E 


Low density 


Gonzalez 


triangles 


NGC 7562 


248. 


.0 


2 


.9 


1 


.69 


0.05 


4.54 


0.04 


2.87 


0.04 


E 


High density 


Gonzalez 


circles 


NGC 7619 


300. 


,3 


3 


,1 


1 


.36 


0.04 


5.06 


0.04 


3.06 


0.04 


E 


High density 


Gonzalez 


circles 


NGC 7626 


253. 


,1 


3 


.0 


1 


.46 


0.05 


5.05 


0.04 


2.83 


0.04 


E 


High density 


Gonzalez 


circles 


NGC 7785 


239. 


,6 


3 


.1 


1 


.63 


0.06 


4.60 


0.04 


2.91 


0.04 


E 


Low density 


Gonzalez 


circles 


NGC 4874 


259. 


,6 


4, 


.3 


2 


.13 


0.02 


4.93 


0.02 


2.91 


0.02 


cD 


High density 


Mehlert 


diamonds 


NGC 4889 


358. 


,7 


4, 


,2 


1 


.93 


0.03 


5.45 


0.04 


3.04 


0.03 


cD 


High density 


Mehlcrt 


diamonds 


NGC 4839 


275. 


,5 


5 


.8 


1 


.42 


0.04 


4.92 


0.04 


2.75 


0.04 


E 


High density 


Mehlert 


circles 


NGC 4841A 


263. 


,9 


7 


.5 


1 


.53 


0.05 


4.51 


0.05 


2.89 


0.04 


E 


High density 


Mehlcrt 


circles 


NGC 4926 


273. 


,3 


4, 


.6 


1 


.50 


0.06 


5.17 


0.06 


2.50 


0.05 


E 


High density 


Mehlcrt 


circles 


IC 4051 


258. 


,7 


2 


,1 


1 


.42 


0.06 


5.34 


0.07 


2.75 


0.05 


E 


High density 


Mehlcrt 


circles 


NGC 4895 


207. 


,2 


1 


,2 


1 


.68 


0.05 


4.39 


0.06 


2.60 


0.05 


SO 


High density 


Mehlcrt 


circles 


NGC 4860 


280. 


.5 


4, 


.1 


1 


.39 


0.06 


5.39 


0.07 


2.85 


0.05 


E 


High density 


Mehlert 


circles 


NGC 4896 


209. 


,1 


4 


.8 


1 


.72 


0.06 


4.04 


0.06 


2.67 


0.05 


SO 


High density 


Mehlert 


circles 


NGC 4865 


249. 


,7 


4, 


,1 


2, 


.26 


0.08 


4.51 


0.09 


3.01 


0.07 


SO 


High density 


Mehlert 


diamonds 


NGC 4923 


186. 


.0 


1, 


.9 


1 


.70 


0.05 


4.43 


0.05 


2.69 


0.04 


E 


High density 


Mehlert 


circles 


GMP 2413 


180. 


.0 


3. 


.0 


1 


.80 


0.06 


4.08 


0.06 


2.60 


0.05 


SO 


High density 


Mehlert 


circles 


NGC 4840 


216. 


.6 


2. 


,7 


1 


.63 


0.07 


4.94 


0.07 


2.91 


0.06 


E 


High density 


Mehlert 


circles 


NGC 4869 


188. 


,1 


2 


.1 


1 


.40 


0.05 


4.83 


0.05 


2.90 


0.04 


E 


High density 


Mehlert 


circles 


NGC 4908 


192. 


,4 


2 


,7 


1 


.58 


0.09 


4.58 


0.09 


2.65 


0.07 


E 


High density 


Mehlcrt 


circles 


IC 4045 


167. 


,3 


2 


,1 


1 


.46 


0.06 


4.70 


0.07 


2.77 


0.05 


E 


High density 


Mehlcrt 


circles 


NGC 4871 


154. 


,3 


3 


.0 


1 


.61 


0.10 


4.36 


0.11 


2.61 


0.09 


SO 


High density 


Mehlcrt 


circles 


NGC 4850 


155. 


,8 


2 


.0 


1 


.57 


0.06 


4.39 


0.06 


2.58 


0.05 


E 


High density 


Mehlert 


circles 


NGC 4883 


157. 


,9 


2 


.0 


1 


.58 


0.06 


4.44 


0.06 


2.86 


0.05 


SO 


High density 


Mehlcrt 


circles 


GMP 1853 


160. 


,4 


5, 


.0 


1 


.63 


0.11 


4.26 


0.12 


2.92 


0.09 


SO 


High density 


Mehlcrt 


circles 


GMP 3661 


154. 


,2 


4, 


.1 


2 


.25 


0.08 


4.20 


0.08 


2.83 


0.07 


SO 


High density 


Mehlcrt 


squares 


GMP 4679 


82, 


,2 


2 


.5 


2 


.39 


0.08 


3.31 


0.09 


2.50 


0.07 


SO 


High density 


Mehlcrt 


triangles 


NGC 4872 


171, 


,7 


3 


.1 


2 


.05 


0.05 


4.05 


0.06 


2.82 


0.04 


E 


High density 


Mehlcrt 


squares 


IC 4041 


114, 


,2 


3, 


.2 


2 


.19 


0.08 


3.82 


0.09 


2.82 


0.07 


SO 


High density 


Mehlert 


squares 


NGC 4957 


208, 


.4 


2, 


.0 


1 


.76 


0.03 


4.53 


0.03 


2.93 


0.02 


E 


High density 


Mehlert 


circles 


NGC 4952 


252, 


,6 


1, 


.7 


1 


.71 


0.03 


4.76 


0.03 


2.69 


0.02 


E 


High density 


Mehlert 


circles 


NGC 4944 


144, 


.8 


1, 


.1 


2 


.13 


0.03 


3.83 


0.03 


2.64 


0.02 


SO 


High density 


Mehlert 


squares 


NGC 4931 


150, 


.8 


1, 


.2 


2 


.10 


0.04 


3.88 


0.05 


2.78 


0.04 


SO 


High density 


Mehlert 


squares 


GMP 1990 


208, 


.9 


2, 


.4 


1 


.40 


0.04 


4.78 


0.04 


2.50 


0.03 


E 


High density 


Mehlert 


circles 


NGC 4827 


243, 


.7 


2, 


.1 


1 


.53 


0.03 


4.89 


0.03 


2.80 


0.02 


E 


High density 


Mehlert 


circles 


NGC 4816 


231, 


.1 


2, 


.4 


1 


.61 


0.02 


4.87 


0.02 


2.74 


0.02 


SO 


High density 


Mehlert 


circles 


NGC 4807 


178, 


.5 


2, 


.0 


1 


.81 


0.06 


4.39 


0.06 


2.78 


0.05 


E 


High density 


Mehlert 


circles 
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Name 


Sigma 


Error 


H/3 


Error 


M 


[gb 


Error 


(Fe> 


Error 


Type 


Environment 


Source 


Class 


ESO 244G045 


185. 


1 


3.7 


1, 


.66 


0.09 


4. 


12 


0.08 


2.86 


0.05 


SO 


Low density 


Beuing 


circles 


ESQ 1070040 


147. 


,0 


5.8 


2. 


.24 


0.25 


3. 


,63 


0.16 


2.97 


0.09 


E 


Low density 


Beuing 


squares 


ESO 1370080 


297. 


3 


7.1 


1, 


.62 


0.18 


5. 


16 


0.15 


3.18 


0.09 


SO 


High density 


Beuing 


diamonds 


ESO 1480170 


134. 


5 


7.5 


2, 


.26 


0.52 


3. 


,49 


0.32 


2.58 


0.20 


E 


Low density 


Beuing 


triangles 


ESO 1850540 


277. 


2 


4.3 


1, 


.57 


0.06 


5. 


11 


0.05 


3.07 


0.03 


E 


High density 


Beuing 


circles 


ESO 2210200 


153. 





5.0 


1, 


.90 


0.17 


4. 


,01 


0.13 


2.83 


0.08 


SO 


Low density 


Beuing 


circles 


ESO 3220600 


193. 


2 


3.8 


1, 


.99 


0.17 


4. 


,45 


0.14 


3.08 


0.08 


SO 


High density 


Beuing 


diamonds 


ESO 4230240 


182. 


5 


5.5 


2. 


.06 


0.37 


3. 


,84 


0.24 


2.97 


0.14 


SO 


Low density 


OUveira 


squares 


ESO 4430240 


287. 





7.8 


1, 


.75 


0.15 


5. 


15 


0.13 


3.13 


0.08 


SO 


High density 


Beuing 


diamonds 


ESO 4450020 


255. 


5 


6.4 


1, 


.69 


0.09 


4. 


,76 


0.08 


3.02 


0.05 


SO 


Low density 


Beuing 


circles 


ESO 4940350 


132. 


1 


2.4 


2, 


.20 


0.16 


3. 


,53 


0.12 


2.74 


0.08 


SO 


Low density 


OUveira 


triangles 


IC 1633 


347. 


7 


7.9 


1, 


,73 


0.09 


5. 


,42 


0.08 


2.98 


0.05 


SO 


High density 


Beuing 


diamonds 


IC 3152 


167. 


,8 


3.7 


1, 


.82 


0.19 


3. 


,66 


0.16 


2.89 


0.10 


SO 


Low density 


Beuing 


circles 


IC 3986 


311. 


5 


7.0 


1, 


.87 


0.15 


5. 


,03 


0.12 


3.12 


0.08 


SO 


High density 


Beuing 


diamonds 


IC 4797 


220. 


6 


6.4 


1, 


.92 


0.26 


4. 


,52 


0.18 


2.75 


0.10 


E 


Low density 


Beuing 


circles 


IC 4931 


289. 





6.0 


1, 


.80 


0.08 


4. 


,76 


0.07 


3.06 


0.04 


SO 


High density 


Beuing 


diamonds 


IC 5105 


308. 


1 


5.4 


1, 


.57 


0.06 


5. 


,26 


0.05 


2.37 


0.03 


SO 


Low density 


OUveira 


circles 


NGC 0312 


254. 


8 


4.4 


1, 


.83 


0.09 


4. 


,56 


0.08 


2.48 


0.05 


E 


Low density 


Beuing 


circles 


NGC 0484 


244. 


5 


4.5 


1, 


.60 


0.08 


4. 


,41 


0.07 


2.95 


0.04 


SO 


Low density 


Beuing 


circles 


NGC 0596 


161. 


8 


1.6 


2, 


.12 


0.05 


3. 


,95 


0.04 


2.81 


0.03 


E 


Low density 


OUveira 


squares 


NGC 0636 


178. 


5 


5.0 


1, 


.86 


0.26 


4. 


,38 


0.17 


2.83 


0.09 


E 


Low density 


OUveira 


circles 


NGC 0731 


162. 


1 


2.6 


2, 


.10 


0.08 


3. 


,92 


0.07 


2.86 


0.04 


SO 


Low density 


Beuing 


squares 


NGC 1052 


202. 


6 


1.7 


1, 


.22 


0.04 


5. 


,53 


0.03 


2.77 


0.02 


E 


Low density 


Beuing 


circles 


NGC 1172 


119. 


3 


1.5 


1, 


.82 


0.11 


3. 


,70 


0.08 


2.52 


0.06 


SO 


Low density 


OUveira 


circles 


NGC 1316 


223. 


1 


2.1 


2, 


.07 


0.03 


3. 


,90 


0.03 


2.84 


0.02 


SO 


High density 


OUveira 


squares 


NGC 1340 


175. 


4 


1.5 


2, 


.23 


0.05 


3. 


,99 


0.04 


3.02 


0.03 


SO 


High density 


OUveira 


squares 


NGC 1375 


75. 


,1 


1.5 


2, 


.35 


0.53 


2. 


,68 


0.33 


2.56 


0.20 


SO 


High density 


Beuing 


triangles 


NGC 1395 


250. 





3.0 


1, 


.62 


0.05 


5. 


,21 


0.04 


2.93 


0.03 


E 


Low density 


OUveira 


circles 


NGC 1400 


256. 


2 


3.0 


1, 


.38 


0.05 


5. 


16 


0.04 


2.93 


0.03 


SO 


Low density 


OUveira 


circles 


NGC 1407 


259. 


7 


3.7 


1, 


.67 


0.07 


4. 


,88 


0.06 


2.85 


0.03 


E 


Low density 


OUveira 


circles 


NGC 1411 


144. 





0.9 


1, 


.74 


0.04 


3. 


,90 


0.04 


2.57 


0.03 


SO 


Low density 


OUveira 


circles 


NGC 1549 


203. 


3 


1.4 


1, 


.79 


0.03 


4. 


,39 


0.03 


2.88 


0.02 


E 


Low density 


OUveira 


circles 


NGC 1553 


207. 


2 


1.1 


1, 


.88 


0.03 


4. 


,44 


0.03 


3.09 


0.02 


SO 


Low density 


OUveira 


circles 


NGC 2434 


180. 


,4 


3.9 


1, 


.87 


0.13 


3. 


,72 


0.10 


2.87 


0.07 


E 


Low density 


Beuing 


circles 


NGC 2717 


210. 





6.5 


1, 


.91 


0.18 


4. 


,69 


0.14 


2.86 


0.09 


SO 


Low density 


OUveira 


circles 


NGC 2891 


100. 


8 


0.6 


2, 


.54 


0.20 


3. 


,07 


0.16 


2.33 


0.11 


SO 


Low density 


Beuing 


triangles 


NGC 2986 


282. 


2 


4.2 


1, 


,48 


0.06 


4. 


,97 


0.05 


2.92 


0.03 


E 


Low density 


OUveira 


circles 


NGC 3078 


268. 


1 


4.2 


1, 


.12 


0.09 


5. 


,20 


0.07 


3.16 


0.04 


E 


Low density 


OUveira 


circles 


NGC 3224 


155. 


,8 


2.3 


2. 


,31 


0.14 


3. 


,91 


0.12 


2.92 


0.08 


E 


High density 


Beuing 


squares 


NGC 3923 


267. 


,9 


3.9 


1. 


,87 


0.08 


5. 


12 


0.07 


3.07 


0.04 


E 


Low density 


Beuing 


diamonds 


NGC 4786 


291. 


,6 


5.8 


1. 


,61 


0.09 


4. 


57 


0.08 


2.87 


0.05 


SO 


Low density 


Beuing 


circles 


NGC 5087 


303. 


,3 


5.1 


1. 


,62 


0.11 


5. 


,39 


0.09 


3.06 


0.05 


SO 


Low density 


Beuing 


circles 


NGC 5516 


313. 


,0 


10.1 


1. 


,46 


0.08 


4. 


,81 


0.06 


3.29 


0.04 


SO 


Low density 


Beuing 


circles 


NGC 5791 


271. 


,8 


9.8 


1. 


,60 


0.19 


5. 


,06 


0.15 


3.30 


0.10 


E 


Low density 


Beuing 


circles 


NGC 5903 


209. 


,2 


3.7 


1. 


,68 


0.10 


4. 


,44 


0.08 


2.90 


0.05 


E 


Low density 


Beuing 


circles 


NGC 6305 


158. 


,9 


4.3 


2. 


,03 


0.13 


4. 


,07 


0.10 


2.80 


0.07 


SO 


Low density 


Beuing 


squares 


NGC 6958 


193. 


,9 


2.1 


1. 


,63 


0.05 


3. 


,89 


0.04 


2.75 


0.03 


SO 


Low density 


Beuing 


circles 


NGC 7041 


238. 


,7 


1.0 


1. 


,88 


0.03 


4. 


,70 


0.03 


2.52 


0.02 


SO 


Low density 


OUveira 


circles 


NGC 7185 


101. 


,2 


1.7 


2. 


,35 


0.16 


3. 


,05 


0.13 


2.63 


0.09 


SO 


Low density 


Beuing 


triangles 


NGC 7365 


115. 


,4 


5.5 


2. 


,20 


0.47 


3. 


,25 


0.30 


2.56 


0.18 


SO 


Low density 


Beuing 


triangles 


NGC 7796 


264. 


1 


5.5 


1. 


,51 


0.08 


5. 


,02 


0.07 


2.78 


0.04 


SO 


Low density 


OUveira 


circles 



Note. — Cols. 3, 5, 7, 9 give 1-a errors. Sources for the data (col. 12): IGonzaled <T99ll . iBeiiing et alJ l200l) . Mchlcrt e^D HOO^), 
C. Mendes de OUveira et al. (in preparation). The last column indicates the classification shown by symbol type in Fig. 
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Name 


Age 


Error 


[Z/H] 


Error 


[a 


:/Fe] 


Error 


NGC 0221 


2, 


,4 


0. 


2 


0. 


,152 


0, 


.030 


-0 


1.025 


0. 


.013 


NGC 0224 


7 


.0 


1 


,7 


0. 


,441 





.048 


0, 


.219 


0. 


.017 


NGC 0315 


6 


.6 


1, 


,4 


0. 


397 


0, 


.038 


0, 


.282 


0, 


.019 


NGC 0507 


8 


,1 


1 


,9 


0. 


261 


0, 


.071 


0, 


.246 


0, 


.040 


NGC 0547 


10, 


,7 


1, 


,7 


0. 


,321 


0, 


.049 


0, 


.311 


0, 


.020 


NGC 0584 


2 


.8 


0, 


,3 


0. 


,478 


0, 


.046 


0, 


.223 


0, 


.014 


NGC 0636 


4, 


,4 


0, 


,6 


0. 


376 


0, 


.028 


0, 


.133 


0, 


.013 


NGC 0720 




,4 


2 


,4 


0. 


485 





.083 


0, 


.353 


0, 


.039 


NGC 0821 


8 


.9 


1, 


,2 


0. 


297 





.036 


0, 


.189 


0, 


.015 


NGC 1453 


9 


,4 


1, 


,6 


0. 


380 





.044 


0, 


.256 


0, 


.018 


NGC 1600 


9 


,7 


1, 


,9 


0. 


420 





.046 


0, 


.263 


0, 


.021 


NGC 1700 


2 


.6 


0, 


,3 


0, 


500 





.049 


0, 


.167 


0, 


.014 


NGC 2300 


7, 


.3 


1, 


.5 


0. 


,427 


0, 


.040 


0, 


.279 


0, 


.017 


NGC 2778 


6, 


,4 


1, 


.9 


0. 


.366 


0, 


.049 


0. 


.268 


0. 


.023 


NGC 3377 


3, 


.6 


0, 


.5 


0. 


,253 


0, 


.035 


0. 


.234 


0. 


.014 


NGC 3379 


10, 


.0 


1 


,1 


0. 


,299 


0, 


.036 


0. 


.259 


0. 


.012 


NGC 3608 


8, 


.0 


1, 


.5 


0. 


.334 


0, 


.042 


0. 


.213 


0. 


.016 


NGC 3818 


6, 


,7 


1, 


.9 


0. 


.422 


0, 


.053 


0. 


.263 


0. 


.025 


NGC 4261 


16 


.3 


1 


.8 


0. 


.275 





.035 


0, 


.255 


0. 


.015 


NGC 4278 


12 


.0 


1 


,3 


0. 


225 


0, 


.038 


0, 


.331 


0. 


.015 


NGC 4374 


12 


.8 


1 


,1 


0. 


207 





.032 


0, 


.259 


0, 


.011 


NGC 4472 


9 


,6 


1, 


,4 


0. 


342 





.046 


0, 


.258 


0, 


.021 


NGC 4478 




,3 


1, 


,2 


0. 


344 


0, 


.049 


0, 


.180 


0, 


.021 


NGC 4489 


2 


.0 


0, 


,3 


0. 


225 





.043 


0, 


.078 


0, 


.021 


NGC 4552 


12 


,4 


1, 


,5 


0. 


356 





.034 


0, 


.277 


0, 


.011 


NGC 4649 


14, 


,1 


1, 


.5 


0. 


362 


0, 


.029 


0, 


.296 


0, 


.012 


NGC 4697 


8, 


.3 


1, 


,4 


0, 


.148 


0, 


.043 


0, 


.155 


0, 


.018 


NGC 5638 


9 


.5 


0, 


.9 


0, 


.274 


0, 


.033 


0, 


.239 


0, 


.015 


NGC 5812 


6 


.5 


1, 


,0 


0. 


.432 


0, 


.029 


0, 


.226 


0, 


.014 


NGC 5813 


16, 


.6 


2, 


,2 


0. 


.059 


0, 


.042 


0, 


.270 


0, 


.018 


NGC 5831 


3 


.0 


0. 


,4 


0. 


.504 


0, 


.043 


0, 


.186 


0, 


.017 


NGC 5846 


14, 


,2 


2, 


2 


0. 


.226 


0, 


.051 


0. 


.270 


0, 


.017 


NGC 6127 


12, 


.6 


1, 


.5 


0. 


.264 


0, 


.044 


0. 


.284 


0. 


.017 


NGC 6702 


1, 


,7 


0, 


,1 


0. 


.686 


0, 


.070 


0. 


.153 


0. 


.017 


NGC 6703 


4, 


.8 


0. 


.8 


0. 


.354 


0, 


.034 


0. 


.180 


0. 


.016 


NGC 7052 


13, 


.2 


2, 


,1 


0. 


.263 


0, 


.050 


0. 


.296 


0. 


.021 


NGC 7454 


3 


.7 


0. 


.6 


0. 


.032 


0, 


.044 


0. 


.117 


0. 


.021 


NGC 7562 


8, 


.6 


1, 


.3 


0. 


.282 





.039 


0, 


.219 


0. 


.014 


NGC 7619 


15 


,4 


1, 


,4 


0. 


291 


0, 


.031 


0, 


.233 


0, 


.013 


NGC 7626 


13 


.9 


1, 


,5 


0. 


251 





.038 


0, 


.302 


0, 


.013 


NGC 7785 


9 


,7 


1, 


,6 


0. 


280 


0, 


.049 


0, 


.209 


0, 


.016 


NGC 4874 


2 


,4 


0, 


,1 


0. 


.719 


0, 


.026 


0, 


.353 


0, 


.008 


NGC 4889 


2 


.8 


0, 


,1 


0. 


769 


0, 


.043 


0, 


.388 


0, 


.020 


NGC 4839 


15 


.8 


1, 


,4 


0. 


159 


0, 


.044 


0, 


.300 


0, 


.027 


NGC 4841A 


12 


,2 


1 


,3 


0. 


.178 


0, 


.049 


0, 


.186 


0, 


.024 


NGC 4926 


14, 


1 


1 


,7 


0. 


.177 


0, 


.059 


0, 


.428 


0, 


.037 


IC 4051 


15 


.0 


1, 


,9 


0. 


276 


0, 


.057 


0, 


.374 


0, 


.032 


NGC 4895 


10 


,4 


1, 


,3 


0. 


115 





.054 


0, 


.266 


0, 


.032 


NGC 4860 


15 


.2 


2, 


,4 


0. 


.315 


0, 


.057 


0, 


.350 


0, 


.038 


NGC 4896 


9, 


,7 


1, 


,4 


0. 


.073 


0, 


.056 


0. 


.170 


0. 


.032 


NGC 4865 


2, 


.0 


0. 


,2 


0. 


.775 


0, 


.114 


0. 


.269 


0. 


.035 


NGC 4923 


9, 


,4 


1, 


,1 


0. 


.180 


0, 


.055 


0. 


.249 


0. 


.037 


GMP 2413 


8, 


.0 


1, 


.5 


0. 


.105 


0, 


.061 


0. 


.214 


0. 


.039 


NGC 4840 


9 


,2 


1, 


.8 


0. 


.371 


0, 


.058 


0. 


.273 


0. 


.030 


NGC 4869 


15 


.4 


1 


,7 


0. 


.190 





.055 


0, 


.236 


0. 


.034 


NGC 4908 


12 


.2 


2, 


,5 


0, 


127 


0, 


.084 


0, 


.278 


0, 


.052 


IG 4045 


14, 


.6 


2, 


,1 


0. 


134 





.075 


0, 


.255 


0, 


.044 


NGC 4871 


11 


.9 


2, 


,7 


0. 


069 


0, 


.103 


0, 


.248 


0, 


.070 


NGC 4850 


13 


,1 


1, 


,6 


0. 


037 


0, 


.068 


0, 


.257 


0, 


.055 


NGC 4883 


11 


.2 


1, 


,9 


0. 


182 


0, 


.078 


0, 


.184 


0, 


.037 


GMP 1853 


10, 


.2 


3, 


,0 


0. 


185 


0, 


.112 


0, 


.130 


0, 


.065 


GMP 3661 


2, 


.1 


0, 


,3 


0, 


,543 


0, 


.107 


0, 


.244 


0, 


.040 


GMP 4679 


2, 


.2 


0, 


,3 


0, 


183 


0, 


.058 


0, 


.150 


0. 


.041 


NGC 4872 


3, 


.5 


0, 


,5 


0, 


,341 


0, 


.064 


0, 


.178 


0, 


.039 


IC 4041 


2, 


.5 


0, 


,4 


0, 


,358 


0, 


.070 


0, 


.147 


0, 


.031 


NGC 4957 


6, 


.5 


0, 


.7 


0, 


,351 


0, 


.022 


0, 


.210 


0, 


.011 


NGC 4952 


8, 


.8 


0. 


.8 


0. 


.279 


0, 


.036 


0. 


.316 


0. 


.021 


NGC 4944 


3, 


.3 


0. 


.3 


0. 


.250 


0, 


.033 


0. 


.195 


0. 


.020 


NGC 4931 


3, 


.3 


0. 


.4 


0. 


.306 


0, 


.042 


0. 


.159 


0. 


.023 


GMP 1990 


18, 


.2 


1, 


.6 


0. 


Oil 


0, 


.042 


0. 


.349 


0. 


.038 


NGC 4827 


12 


.0 


0, 


.9 


0. 


.250 


0, 


.041 


0. 


.287 


0. 


.025 


NGC 4816 


10 


.6 


0, 


.5 


0. 


.267 





.030 


0. 


.311 


0. 


.020 


NGC 4807 


6 


.4 


1, 


.5 


0. 


.274 


0, 


.081 


0. 


.232 


0. 


.048 
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Table A3. - continued. 



Name 


Age 


Error 


[Z/H] 


Error 


[q 


:/Fc] 


Error 


ESO 244G045 


9 


.9 


2, 


,1 


0. 


141 


0. 


064 


0, 


.121 


0, 


.025 


ESQ 1070040 


2 


,1 


2, 


.3 


0. 


414 


0. 


190 


0, 


.081 


0, 


.059 


ESO 1370080 


7, 


.0 


4, 


,0 


0. 


517 


0. 


130 


0, 


.251 


0, 


.042 


ESO 1480170 


2 


,7 


4, 


,0 


0. 


196 


0. 


189 


0, 


.153 


0, 


.088 


ESO 1850540 


9 


,2 


1, 


.5 


0. 


434 


0. 


035 


0, 


.259 


0, 


.014 


ESO 2210200 


5 


,1 


2, 


.5 


0. 


246 


0. 


112 


0, 


.148 


0, 


.045 


ESO 3220600 


3 


.0 


1, 


,7 


0. 


533 


0. 


154 


0, 


.192 


0, 


.044 


ESO 4230240 


3 


.3 


4, 


,7 


0. 


324 


0. 


220 


0, 


.125 


0, 


.083 


ESO 4430240 


4, 


,4 


2, 


,4 


0. 


575 


0. 


107 


0, 


.277 


0, 


.040 


ESO 4450020 


7, 


,2 


2, 


,0 


0. 


402 


0. 


060 


0, 


.223 


0, 


.026 


ESO 4940350 


2 


.8 


1, 


,2 


0. 


226 


0. 


097 


0, 


.150 


0, 


.043 


IC 1633 


5 


,2 


1, 


.6 


0. 


563 


0. 


057 


0, 


.366 


0, 


.024 


IC 3152 


6 


,4 


3, 


,6 


0. 


132 


0. 


131 


0, 


.034 


0, 


.051 


IC 3986 


3 


.3 


2, 


,0 


0. 


623 


0. 


131 


0, 


.280 


0, 


.039 


IC 4797 


4, 


,7 


4, 


,2 


0. 


362 


0. 


181 


0, 


.279 


0, 


.067 


IC 4931 


4 


,7 


1, 


,5 


0. 


475 


0. 


057 


0, 


.230 


0, 


.023 


IC 5105 


11 


,4 


1, 


,7 


0. 


210 


0. 


041 


0, 


.481 


0, 


.013 


NGC 0312 


7, 


.3 


2, 


.0 


0. 


205 


0. 


063 


0, 


.362 


0, 


.026 


NGC 0484 


10, 


,4 


2, 


.0 


0. 


222 


0. 


059 


0, 


.155 


0, 


.022 


NGC 0596 


2 


.8 


0, 


,4 


0. 


360 


0. 


035 


0, 


.175 


0, 


.015 


NGC 0636 


4, 


,7 


4, 


,1 


0. 


360 


0. 


188 


0, 


.225 


0, 


.056 


NGC 0731 


3 


,1 


0, 


,7 


0. 


351 


0. 


066 


0, 


.147 


0, 


.026 


NGC 1052 


21 


,7 


1, 


,6 


0. 


222 


0. 


025 


0, 


.390 


0, 


.007 


NGC 1172 


7, 


.0 


2, 


,3 


0. 


077 


0. 


057 


0, 


.088 


0, 


.031 


NGC 1316 


3 


,2 


0, 


,2 


0. 


338 


0. 


023 


0, 


.150 


0, 


.009 


NGC 1340 


2 


,1 


0, 


,1 


0. 


490 


0. 


056 


0, 


.200 


0, 


.014 


NGC 1375 


2 


.6 


3, 


,1 


0. 


Oil 


0. 


177 


-0 


.044 


0, 


.133 


NGC 1395 


7, 


.6 


1, 


,4 


0. 


439 


0. 


033 


0, 


.353 


0, 


.013 


NGC 1400 


15 


.0 


1, 


.6 


0. 


226 


0. 


032 


0, 


.338 


0, 


.010 


NGC 1407 


7, 


,4 


1, 


,8 


0. 


379 


0. 


044 


0, 


.300 


0, 


.019 


NGC 1411 


9 


1 


0, 


,9 


0. 


042 


0. 


026 


0, 


.163 


0, 


.015 


NGC 1549 


6 


,1 


0, 


.6 


0. 


283 


0. 


020 


0, 


.237 


0, 


.009 


NGC 1553 


4, 


,7 


0, 


,3 


0. 


370 


0. 


017 


0, 


.188 


0, 


.009 


NGC 2434 


5 


.5 


2, 


,4 


0. 


167 


0. 


090 


0, 


.066 


0, 


.038 


NGC 2717 


3 


,7 


2, 


.8 


0. 


460 


0. 


128 


0, 


.300 


0, 


.042 


NGC 2891 


2 


.0 


0, 


.8 


0. 


113 


0. 


099 


0, 


.163 


0, 


.057 


NGC 2986 


11 


,2 


1, 


,9 


0. 


347 


0. 


045 


0, 


.250 


0, 


.014 


NGC 3078 


22 


.3 


4, 


,2 


0. 


228 


0. 


055 


0, 


.250 


0, 


.016 


NGC 3224 


1 


.9 


0, 


,7 


0. 


554 


0. 


139 


0, 


.166 


0, 


.041 


NGC 3923 


3, 


.3 


0, 


,8 


0, 


,623 


0, 


,070 


0, 


.308 


0, 


.022 


NGC 4786 


10, 


.4 


2, 


,4 


0, 


,237 


0, 


,063 


0, 


.213 


0, 


.024 


NGC 5087 


7, 


.3 


2, 


.3 


0, 


,523 


0, 


,068 


0, 


.323 


0, 


.026 


NGC 5516 


11, 


.2 


2, 


,1 


0, 


,380 


0. 


,046 


0, 


.131 


0, 


.017 


NGC 5791 


6, 


.8 


3, 


,9 


0, 


,530 


0. 


124 


0, 


.199 


0, 


.045 


NGC 5903 


8, 


.8 


2, 


.4 


0, 


,264 


0, 


,075 


0, 


.186 


0, 


.029 


NGC 6305 


3, 


.7 


1, 


.6 


0, 


,322 


0, 


,092 


0, 


.188 


0, 


.035 


NGC 6958 


11, 


.3 


1, 


,0 


0, 


,015 


0, 


,027 


0, 


.097 


0, 


.013 


NGC 7041 


4, 


.6 


0, 


.3 


0, 


,373 


0, 


,025 


0, 


.350 


0, 


.010 


NGC 7185 


2, 


.3 


0, 


.9 


0, 


145 


0, 


100 


0, 


.038 


0, 


.047 


NGC 7365 


3, 


.1 


3, 


,4 


0, 


,090 


0, 


176 


0, 


.095 


0, 


.104 


NGC 7796 


11, 


.8 


1, 


,9 


0, 


,248 


0, 


,050 


0, 


.344 


0, 


.019 



Not e. — The stellar population parameters age, total mctallicity [Z/H], and [a/Fo] abundance ratio are derived from H/3, Mgb, and (Fe) 
(Table ESJ using SSP models of Thomas ct ah 1 2003 b) and Th omas ct ah ( 2004). Cols. 3, 5, 7 give l-cr errors. The latter are estimated 
via Monte Carlo simulations (see text). 
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Mg b Mg b 

Fig. Bll. — Lick indices [MgFe]' versus H/3 (top panels) and Mgb versus (Fe) (bottom panels). Symbols like in Fig.^l SSP models with 
the metallicities [Z/H] = 0.0, 0.35, 0.67 are plotted for different ages i = 2, 3, 5, 10, 15 Gyr at fixed a/Fe ratio ([a/Fe] = 0) (top panels) and 
different a/Fe ratios [a/Fe] = 0.0, 0.3, 0.5 at fixed age (t = 12 Gyr) (bottom panels) as indicated by the labels. Different from Fig. l2l 
the models are corrected for the decrease of O/Fe with increasing metallicity at super-solar metallicity in the solar neighborhood assuming 
oxygen and the other a-elements being lumped together. 



THE LOCAL ABUNDANCE PATTERN AT HIGH METALLICITIES 

In a recent study. IProctor et al.l l)2004() argue that the correlation of a/Fe with velocity dispersion is an 'artefact 
of incomplete calibrat ion (of the models) to the Lick system'. The basis for their conclusion is the recent finding 
of iBensbv et al.l l)2004j) . that the Milky Way disk stars with super-solar metallicities do not have solar but sub-solar 
O/Fe ratios, O/Fe decreasing with increasing Fe/H (see also lAllende Prieto et al.ll2004l: IPompeia et al.ll2003l ). Stellar 
population models, being calibrated with these stars, are then likely to reflect sub-solar O/Fe ratios at high metallicities. 
If the element Oxygen is taken as representative for all a-elements, this would obviously lead to an overest imation of the 
O/Fe ratio (or a/Fe ratio) in stellar populations with super-solar metallicities like early- type galaxies. IProctor et alJ 
l)2004j) show that the correction for this effect makes the a/Fe-sigma relation disappear. 

However, there is a further very important (and somewhat puzzling) point shown in IBensbv et al.l l)2004() . All 
the other a-elements, Mg, Si, Al, Na, Ca, and Ti, do not sh ow this pattern, but level-off at a solar X/Fe ratio in 
stars wi th super-solar metallicities jB cnsbv ct al. 200^1 iFuhrma nn 1998). Only the element oxygen (and maybe also 
Carbon, iThoren. Edvardsson. & Custafsson .20041 sticks out of this homogeneous behavior. This differential pattern 
between Oxygen and the other a-elements cannot be neglected in models taking the local abundance pattern at high 
metallicities into account. This is particularly important in the context of this work, as the a/Fe ratios presented here 
are der ived from the absorption line indices Mg6 and (Fe), which are sensitive mostly to the elements Fe and Mg, but 
not O l|Korn et a,].ll2nn,^ ITHnicco fc''Be1]ll1 99,^ . 

We investigated this issue constructing two additional flavors of stellar population models that account for the local 
abundance pattern at high metallicity. One in which Oxygen is tak en as representative for all a-elements assuming 
these elements to be enriched in lockstep (like in IProctor et al.ll2004D . and one in which Oxygen is detached from the 
rest of the a-elements as discussed ab ove. The decrease o f O/Fe with metallicity is adopted directly from .Bcnsbv ct all 
l)2004ll and can be parameterized as (jProctor et al.ll200^ : 



[O/Fe] = -0.5 [Fe/H] for [Fe/H] > 



(Bl) 
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Fig. B12. — Stellar population parameters as a functions of velocity dispersion (measured within 1/10 r^) a nd ste llar mass (upper x-axis) 
for the 'old' subpopulation (see Fig.|HJ. Ages and element abundances are derived with the SSP model of Fig. lBlfl Black points are Monte 
Carlo simulations asuming Egn. [Tl (solid lines), which reproduce the data in Fig. 151 



Oxygen and a-elements in lockstep 

The resulting model in the index-index planes, assuming oxygen as a representative for all a-elements, is shown in 
Fig. iBlll (equivalent to Fig. 01. 

The effect is striking. As oxygen dominates metallicity, the correction for the sub-solar 0/Fe of the calibrating stars 
at high metallicities has the effect that the models are practically 'labeled' with lower metallicities. As a consequence, 
stronger metal indices at a given metallicity a re pr edicted, which leads to the derivation of somewhat lower (but 
still super-solar) metallicities (top panel of Fig. IBll'jl . The derived ages, instead, remain the unchanged. Even more 
dramatic is the modification of the model in the Mg6-(Fe) plane (bottom panel). At super-solar metallicity, the 0/Fe 
(hence a/Fe in this case) of the original model based on the local abundance pattern has to be increased in order to 
obtain solar-scaled abundance ratios. This is achieved mainly by decreasing the abundance of iron, so that the iron 
indices get weaker and the Mg h index get stronger (due to the inverse response to Fe abundance, see TMB for details). 
The model slope in this diagram is now close to the slope of the data, which clearly impacts on the derived a/Fe-u 
relation. 

For a better illustratio n, the stellar population parameters of the galaxy sample investigated here obtained with this 
model are plotted in Fig. lB12l fthe equivalent to Fig.|HJl. The relationships of Eqn.Qplus the Monte Carlo realizations 
of Fig. El are also shown for comparison. The a/Fe ratios are lower by about 0.1 dex, and in particular the a/Ye-a 
relation disappears. Metallicities are also lower by about 0.1 dcx, and the Z/H-cr relationship gets significantly tighter. 
These findings are in agreement with the conclusions of Proctor et al. (2004). Note that the ages do not change. 

Oxygen and a-elements separated 

In this second flavor of the model we account for the fact that the other a-elements (in particular magnesium) 
actually do not foll ow the decrease of the local 0/Fe ratio with increasing metallicity at super-solar metallicities 
ijBensbv et al.ll2nn4D . hence [Ne,Mg,Si,S,Ar,Ca,Ti/Fe]= for [Fe/H] > 0. The impact on the index-index diagrams is 
shown in Fig. IB 131 The model is almost indistinguishable from the original model used for the analysis of this pape r 
(Fig.|2Il and very different from the model in which oxygen and the other a-elements are lumped together (Fig. IBlTll . 
The reason is the following. Oxygen contributes almost half, and the rest of the a-elements about one quarter to 
total metallici ty. A variation of the a/Fe ratio including oxygen is therefore obtained mainly by a modification of Fe 
abundance (,Trager et al.ll2000bl and TMB). 

The correction from the local sub-solar 0/Fe to solar element ratios is therefore achieved by a decrease of Fe 
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Mg b Mg b 

Fig. B13. — Lick indices [MgFe]' versus H/3 (top panels) and Mgb versus (Fe) (bottom panels). Symbols like in Fig.^l SSP models with 
the metallicities [Z/H] = 0.0, 0.35, 0.67 are plotted for different ages i = 2, 3, 5, 10, 15 Gyr at fixed a/Fe ratio ([a/Fe] = 0) (top panels) and 
different a/Fe ratios [a/Fe] = 0.0, 0.3, 0.5 at fixed age (t = 12 Gyr) (bottom panels) as indicated by the labels. Different from Fig. l2l 
the models are corrected for the decrease of O/Fe with increasing metallicity at super-solar metallicity in the solar neighborhood assuming 
oxygen and the other a-elements being separated. 



abundance by 40 per cent. When oxygen is detached from the other a-elements, instead, the picture is very different. 
The required decrease of Fe abundance is only 20 per cent, which reduces both the decrease of the Fe indices and 
the increase of Mg6 (note that Mg6 responses inversely to Fe abundance). On top of this, the abundances of the 
other a-elements (including magnesium) need to be reduced by the same amount to maintain the solar a/Fe ratio. 
This leads t o a further decrease of Mg b and an incr ease of the Fe indices, the latter being inversely coupled with Mg 
abundance l|Korn et al.ll2005t iTripicco fc BeU"1995!). The effect is sufficient to bring the model back to its original 
position. As a consequence, the stellar populati on pa rameters derived with this model are excellent agreement with 
the ones presented in this work as shown in Fig. IB 141 

To conclude, when the non-solar 0/Mg ratios of metal-rich stars in the solar neighborhood are taken into account, 
the correction for the local abundance pattern has no significant impact on the model. This implies that the results 
of this paper are robust against this problem. In particular, the a/Fe-a relation found here is real and not caused by 
calibration artefacts. 
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Fig. B14. — Stellar population parameters as a functions of velocity dispersion (measured within 1/10 r^) a nd ste llar mass (upper x-axis) 
for the 'old' subpopulation (see Fig.|HJ. Ages and element abundances are derived with the SSP model of Fig. lB13l Black points are Monte 
Carlo simulations asuming Egn. [Tl (solid lines), which reproduce the data in Fig. 151 



